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Zusammenfassung 



Kontext 

Wochentlich werden neue Entdeckungen von Planeten publiziert, die mittels Transitme- 
thode gewonnen wurden. Irdische Beobachtungen sowie Weltraummissionen wie CoRoT 
und Kepler sind dafur verantwortlich, dass wir die statistischen Liicken mit Planeten 
von entfernten Sternensystemen fiillen. Das ultimative Ziel ist die Entdeckung eines 
habitablen Planeten, vielleicht einer zweiten Erde. 



Ziele 

Ich mochte die stellare Aktivitat und ihren Einfluss auf die Entdeckung von extrasolaren 
Planeten erortern. Bis heute bilden CoRoT-7b und Kepler- 10b die beiden Ausnahmen 
als kleine Gesteinsplaneten genannt „ Super- Erden". Die Frage wirft sich auf, warum un- 
ter den iiber 500 entdeckten und verifizierten Planeten die Anzahl der kleinen Planeten 
derart gering ist. Unsere Statistik die wir bis heute iiber die Verteilung der planeta- 
ren Massen und Durchmesser gesammelt haben, ist durch die Beobachtung verzerrt. 
Ein anderer Grund dafur konnte, abgesehen von der Schwierigkeit kleine Planeten zu 
verifizieren, der hohe Level an stellarer Aktivitat sein, der bisher beobachtet wurde. 

Stellare Aktivitat verlauft auf unterschiedlichen Zeitskalen, von langjahrigen Strah- 
lungsanderungen wie dem bekannten Sonnenzyklus, iiber die stellare Rotation im Be- 
reich von mehreren Tagen bis hin zur Beobachtung von akustischen Moden im Bereich 
von Minuten. Aber auch aperiodische Vorgange wie Flares oder das Aktivitatssignal 
der Granulation kann die Entdeckung eines Planetentransits verhindern. 
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Methoden 

Ich mochte unterschiedliche Methoden zur Detektion von transit-ahnlichen Signalen 
beschreiben, die unterschiedlichen Einfliisse des Sterns, des beobachtenden Instruments 
und deren Auswirkung auf den Erfolg der Methodik. Unterschiedliche Filter-Methoden 
werden diskutiert um die Moglichkeit einer Transitdetektion zu vergroBern und um das 
Transitsignal vom intrinsischen stellaren Signal zu trennen. 

Letztendlich werden unterschiedliche mathematische Modelle und Approximationen 
von Transitfunktionen auf deren Empfindlichkeit auf stellare Aktivitat hin untersucht. 

Ergebnisse 

Ein statistischer Uberblick iiber die stellare Aktivitat im CoRoT-Datenarchiv wird pra- 
sentiert. Der Einfluss stellarer Aktivitat auf unterschiedliche Transitplaneten wird bei 
CoRoT-2b, CoRoT-4b und CoRoT-6b untersucht. 

Schlussfolgerung 

Stellare Aktivitat kann die erfolgreiche Entdeckung eines Planetentransits verhindern, 
wobei CoRoT-7b sicherlich die Grenze markiert. Zukiinftige Missionen wie Plato sind 
notwendig um langjahrige Beobachtungen mit einer Genauigkeit im mmag Bereich zu 
liefern um die Einschrankungen zu iiberwinden, die durch die aktiven Sterne in unserer 
galaktischen Nachbarschaft gegeben sind. 
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Summary 



Context 

New discoveries of transiting extrasolar planets are reported weekly. Ground based 
surveys as well as space borne observatories like CoRoT and Kepler are responsible for 
filling the statistical voids of planets on distant stellar systems. The ultimate goal is 
the discovery of a habitable planet and maybe a second Earth. 

Aims 

I want to discuss the stellar activity and its impact on the discovery of extrasolar planets. 
Up to now the discovery of small rocky planets called "Super- Earths" like CoRoT-7b and 
Kepler-lOb are the only exceptions. The question arises, why among over 500 detected 
and verified planets the amount of smaller planets is strikingly low. Our statistics, we 
collected so far on the distribution of planetary masses and radii, is obviously biased by 
observational effects. Another explanation besides that the verification of small planets 
is an intriguing task, is the high level of stellar activity that has been observed. 

Stellar activity can be observed at different time-scales from long term irradiance 
variations similar to the well known solar cycle, over stellar rotation in the regime of 
days, down to the observations of acoustic modes in the domain of minutes. But also 
non periodic events like flares or the activity signal of the granulation can prevent the 
detection of a transiting Earth sized planet. 

Methods 

I will describe different methods to detect transit-like signals in stellar photometric 
data, the different influences introduced by the star, the observer and their impact 
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on the success. Different filtering techniques will be discussed to improve the transit 
detection capability and to separate the transit signal from the intrinsic stellar signal. 

Finally different mathematical models and approximations of transit signals will be 
examined on their sensibility of stellar activity. 

Results 

We present a statistical overview of stellar activity in the CoRoT dataset. The influence 
of stellar activity will be analysed on different transiting planets: CoRoT-2b, CoRoT-4b 
und CoRoT-6b. 

Conclusions 

Stellar activity can prevent the successful detection of a transiting planet, where CoRoT- 
7b marks the borderline. Future missions like Plato will be required to provide long-term 
observations with mmag precision to overcome the limitations set by active stars in our 
Galactic neighbourhood. 
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1. Introduction 



With the invention of the telescope the last unknown planets of our solar system were 
detected by Herschel in 1781, who discovered Uranus, and Le Verrier (1845), who 
predicted the position of Neptune based on observations of Uranus. It took again 
about 150 years for the detection of the first extrasolar planet (or "exoplanet" for short) 
around the main sequence star 51 Pegasi by Mayor et al. (1995). The first detection of 



a transiting planet came shortly after by Charbonneau et al. (2000). 

Since then the number of detected exoplanets has been increasing rapidly over the 



past years. This fact was also supported by dedicated space missions like MOST (Rowe 



et al, 2006), CoRoT (Baglin et al. 2009) and Kepler (Borucki, 2011) as well as ground 



based surveys like the HAT- network (Bakos et al. , 2002) or WASP and Super- WASP 



(Pollacco et al. , 2006). The technological progress has also improved the detection 



capabilities of radial velocity searches by e.g. the HARPS spectrograph (Pepe et al. 



2004) and SOPHIE (e.g. Bouchy et al. 2009). The large observatories VLT and the 



upcoming E-ELT will show us the next step of direct imaging (Kalas et al. , 2009). 

One the long path to direct observation we are still tied to transit- and radial- velocity- 
measurements of the extrasolar planets to derive the important parameters like mass 



and radius. Upcoming missions like Plato (Catala, 2009) will observe with unpreced- 
ented precision, revealing even more details, but also dealing with new problems. Even 
then the stellar signal becomes more and more relevant. The characterization and the 
understanding of the stellar activity on planet hosting stars will help to explain the 
planetary evolution and its habitability. 

Spectro-polarimetric observations of transits will be required to answer the questions 
of habitability. The interaction of the planets' atmosphere with the stellar wind of the 
host star can only be observed in wavelengths of the respective atmospheric species. 

By combination of the transit method with the observation in different wavelengths, 
the transiting planet with its opaque disk "samples" the surface of the star. This will 
reveal the host star's activity in its photosphere, chromosphere and even corona. If 
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we imagine an ultra-high precision photometry of a small transiting body, we can even 
think of observing the granulation or corona of a distant star. 

So far no working definition for an extrasolar planet has been agreed upon. The 
latest definition of planets dates back to 2006 in the IAU Resolution 5A 2006^j 

RESOLUTION 5A 

The IAU therefore resolves that planets and other bodies in our Solar Sys- 
tem, except satellites, be defined into three distinct categories in the follow- 
ing way: 

1. A planet is a celestial body that 

a) is in orbit around the Sun, 

b) has sufficient mass for its self-gravity to overcome rigid body forces 
so that it assumes a hydrostatic equilibrium (nearly round) shape, 
and 

c) has cleared the neighbourhood around its orbit. 

2. A dwarf planet is a celestial body that 

a) is in orbit around the Sun, 

b) has sufficient mass for its self-gravity to overcome rigid body forces 
so that it assumes a hydrostatic equilibrium (nearly round) shape, 

c) has not cleared the neighbourhood around its orbit, and 

d) is not a satellite. 

e) All other objects , except satellites, orbiting the Sun shall be re- 
ferred to collectively as Small Solar-System Bodies. 

This resolution does not apply to extrasolar planets and has been refined in the 



working definition 2006 (Boss et al. 2007): 
[A]n exoplanet 

• is an object with a true mass below the limiting mass for thermonuclear 
fusion of deuterium (currently calculated to be 13 Jupiter masses for 
objects of solar metallicity), 

• is in orbit around a star or stellar remnant, 



J http : //www. iau. org/ static/archives/releases/pdf /iau0603 . pdf 
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has a mass and/or size that is superior to the one used as a limit for a 
planet in our Solar System. 



We see that the limitations regarding the mass and the size are rather weakly defined. 
Especially the limiting mass for thermonuclear fusion of deuterium is estimated at 
~ 13.0 ± 0.8Mj up , but depends on the pre-conditions like metallicity and can vary 
from « 11.0Mj U p to 16. 3Mj up ( Spiegel et al. 2011). The definition also includes objects 



orbiting Pulsars, but excludes free floating planets (Han, 2006, and references therein). 



1.1. Stellar and Solar Activity 



Stellar activity in the context of transiting exoplanets has been investigated by Alapini 



and Aigrain (2009) for the radial velocity method (see Section 1.2.2). 



Stellar activity happens on different scales and has many varieties. We are interested 
in activity signals that have an impact in extrasolar planet detection and have an 
implication on the evolution and the environment of those planets. 

At first we take a look at the solar activity that has been studied for centuries since 
the detection of the first sunspots. Once again the sun is an ideal template - and in 
frame of the work - for the stars hosting exoplanets. Attempts to deduce stellar activity 



levels from solar variations have been carried out by e. g. Carpano and Fridlund (2008 ). 



The main drawback is the fact that the solar activity levels only correspond to a stellar 
analogue that matches the spectral type and the age of the sun. So far transit search 



programs only cover a fraction of a solar/stellar cycle. The star HD 49933 (Garcia et al. 



2010) shows evidence for magnetic activity that might vary over longer time scales. 



1.1.1. Solar Activity 

The solar activity has been known since the invention of the telescope 400 years ago. 
Spots on the sun have been observed even longer. The solar rotation as observed by the 
virtual motion of the solar spots across the solar surface is a signature of the magnetic 
activity of our central star. As the occurrence varies over the solar cycle, other non 
periodic events like flares occur along with it. 

The solar oscillation namely the p-modes were discovered on the sun, revealing a new 
method to actually probe the interior of our sun. An overview of the solar activity is 
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shown in Figure 1.1 



V 




solar activity 




Figure 1.1.: Mindmap of solar activity. 



Virgo/SoHo-data has been used in the past to simulate realistic lightcurves for transit 



based observations by Aigrain et al. (2003), where different time-scales r of solar activ- 



ity were found (see Table 1.1) 



The power-spectrum P{v) of the solar flux can be described by power-laws for non- 
periodic 



PM 



1 + (2vrz/r) 6 



'1.11 
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1.1. Stellar and Solar Activity 



Table 1.1.: Activity sources and their respective timescales r taken from Aigrain et ah 
(2003) and frequencies v. 



activity source 



timescale r frequency v 



active regions 
sup er-gr anulat ion 
meso-gr anulat ion 
granulation 
bright points 



1...3 x 10 5 s 
3... 7 x 10 4 s 
8000 s 
200 ...500s 
70s 



3 . . . 10 M-Hz 
14 ... 33 M-Hz 
125 M-Hz 
2 ... 5 mHz 
14mHz 



and periodic components 



PM 




[1.2) 



where A is the amplitude, b is the decay time and c represents the frequency trend for 



a central frequency z/ (Harvey et al. 1993). 



The known solar activity signals can be separated in two large groups, periodic and 
non-periodic or eruptive solar activity. 



p-mode Oscillations 

A small scale periodic variability are the 5-minutes oscillations, which were later called 
solar p-modes and were firstly observed in wide band photometry with the ACRIM 



instrument (Woodard and Hudson, 1983). The explanation for the observation is the 



change in temperature of the photosphere (Renaud et al. , 1999). The discovery of these 



acoustic modes goes back to 1960 (Stix, 2002), where periodic changes were visible 
in "Doppler plates". Their maximum amplitude is in the regime of 2 ... 5 mHz which 
roughly corresponds to 3. . . 8 minutes. 



Solar Irradiance 



The Total Solar Irradiance (TSI) has been measured with a high precision, 10 ppm 



according to Willson and Mordvinov (1999) over two cycles. We need continuous space 



based observations to compare the solar activity with the stellar activity on the corres- 
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ponding time scales. The variation with the solar cycles can be seen in Figure L2fl 



TOTAL SOLAR IRRADIANCE MONITORING RESULTS: 1978 to Present 
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Figure 1.2.: Total Solar Irradiance measured with different instruments. The variation 
of the activity levels during the solar cycle is clearly visible. The offsets 
between the different missions arise from unknown instrumental origin. 

TSI data is taken from the SORCE instrument^ The total irradiance is modulated 
by the solar rotation, being visible due to the active regions virtually moving across the 
solar disk and the magnetic cycle, that triggers the appearance of solar activity. 



Solar Rotation and Sunspots 



As mentioned in Willson and Mordvinov (1999) the spectral power at 27 days is evident 



which corresponds to the solar rotation period. The inclination angle % was initially 



2 Data taken from http://www.acrim.com/ 



3 http : / /lasp . Colorado . edu/ sorce/data/tsi_data. htm 
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defined by Carrington (1863) [sic] and measured at Kanzelhohe, Austria i = 7.137° ± 



0.017° (Stix, 2002). 



Solar Flares and Active Regions 

Solar flares also contribute to the TSI. Only the most energetic flares are visible in the 
continuum, otherwise the most energy in the visible spectrum is emitted in Ha. The 
most prominent example of a white light flare was the Haloween event on 28 Oct. 2003 



(Kopp et al., 2004 Woods and Kopp, 2005). 



1.1.2. Stellar Activity 

Taking the sun as a template for other stars hosting planets, we might expect to observe 



the same phenomena as on the sun we described in Section |1.1.1| The main difference 
is the big variety, that we have to expect in the observation of numerous targets. Planet 
hosting stars will have different spectral types than the sun, but also different evolu- 



tionary stages and by that different levels of activity (Guinan et al. 2003). The size 



of the star as listed in Table 1.2, has a big influence on the detected transit depth, 



which will be discussed in Section 1.3 The inclination axis of the host star plays an 
important role, but especially for transit observations, where the observed inclination 
of the system is close to 90°, the expected inclination of the host star's rotation axis 
can be assumed to have a similar value. 

Stellar activity is the limiting factor in the detection and verification of planetary 



candidates (Dumusque et al., 2011). If we want to derive stellar activity levels from the 



sun we need to scale at least with the stellar radii (Aigrain et al. , 2004) 



Stellar activity can be characterized in the power spectrum of the light curve as a 
'generalized' lorentzian profile 



PM 



N 

§1 + 



A 



+ B, 



;i.3) 



where P(v) is the spectral power, on = 2 if a temporal phenomenon is observed, A is the 



amplitude and 1/Az/j is the characteristic timescale (Hulot et al. 2011, and references 



therein). Equation 1.3 is closely related to Equation 1.2 
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Table 1.2.: Spectral types of dwarf stars taken from Cox (2000) with their respective 
masses, radii and effective temperatures. 
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Oscillations 

Solar like oscillations in stars are difficult to detect photometrically, they have been at 



least proven by the radial velocity measurements (Teixeira et al. , 2008; Eggenberger 



et al. , 2004). The detection of p-modes on Procyon is controversial (Matthews et al. 



2004; Regulo and Roca Cortes, 2005). Stellar oscillations have been measured by differ- 



ent space-based missions e.g. CoRoT (Michel et al. 2008) and Kepler (Gilliland et al. 



2010) 



Even observations of p-modes have been reported on planet hosting stars (Ballot 



et al. 2011 Gilliland et al., 2010). Low-degree p-modes were successfully observed by 



CoRoT on the star HD 49933, which is of F5V type for 60 days(Appourchaux et al. 



2008). The maximum amplitude of 3.75 ± 0.23 ppm was found at 1760 \xtlz. 



As reported by Baudin et al. (2011 ) the average amplitude of the oscillations is around 
75 ppm with maximum amplitudes up to several 100 ppm. The maximum amplitude is 
correlated with z/ max , where the lowest frequencies about 20 \jlHz corresponds to amp- 
litudes of 150 ppm. The highest frequency at z/ max w 60 \xH.z has the lowest amplitude 
of 40 ppm. 
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The F5V star HD 49933 shows solar like oscillations that have been measured by 



CoRoT (Appourchaux et al. , 2008). Also the Kepler mission was able to observe acous- 



tic modes (Mathur et al. 2011). The amplitudes can reach 200 ppm in a range of 20 



to 60 |j.Hz(Baudin et al., 2011). 



Stellar Rotation 

Stellar rotation is similar to the solar rotation. It can be identified due to the movement 
of active areas across the stellar disk modulating the total flux that can be measured. 
But different inclination angles have to be taken into account. The inclination angles 
are uniformly distributed, so we also have to expect top-on observations. 




20 



40 



60 80 
time [d] 



100 120 



140 



Figure 1.3.: Stellar rotation as a sinusodial modulation of the flux. The highest amp- 
litude is shown in green. 



Flare Stars 

About 1101 active dwarf stars of spectral type F, G and K are known at the time 
being (July 201 1)Q Among them only 57 flare stars are known. For example the star 

4 http : / / simbad . u- strasbg . f r/ simbad/sim- sam 
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V* VI 041 Tau belongs to the same spectral class G2V as the sun and is a member of 



the Pleiades field. According to Gorlova et al. (2006) this star has a spectral type of 



G9. Nevertheless this object seems to be an interesting target. The star's properties 
are listed in table II .31 

Table 1.3.: Properties of V* V1041 Tau as an example for a solar like Hare star. 



Property Value 
Identifier V* V1014 Tau 

HII 0738 

2MASS J03453940+2345154 
Spectral Type G2V (G9) 

B 13.42 
V 12.26 
J 9.767 
H 9.224 
K 9.007 
(V-K)o 2.30 



So far only little evident is known for flare activity at transiting planets. One of 



the few exceptions is a flare that has been observed at OGLE-TRlOb by |Bentley et al. 
(120091. 



Blend Scenarios 

Since binary stars are very common throughout the galactic neighbourhood, we must 
expect them in transit searches. More or less they are closely related to transiting 
exoplanet systems, where only the secondary body is of smaller size. Brown dwarfs, 
however fall into this regime. The binarity is no intrinsic stellar activity, but definitely 
obfuscating the transit detection. 

There are many different cases that can mimic a transiting planet. Most of them 



are caused by binary stars (e.g. Moutou et al. 2006). These cases are illustrated in 



Figure |1.4| The triple system does not have to be physically connected. Usually a 
background binary contaminates the PSF. These systems can be identified with high 
resolution imaging or spectroscopy. 
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Figure 1.4.: Blend scenarios adopted from (Moutou et al., 2006). a) Grazing binary; 

b) small stellar companion; c) eclipsing binary in a triple system; false 
positive caused by stellar activity or instrumental noise. 



1.2. Detection of Exoplanets 

Many different ways of detecting extrasolar planets have been developed through the 
years, among them the radial-velocity or Doppler method that initially came from 
the investigation of close binaries. The transit method that evolved on the ground 
with several surveys finally went to space leading to successful mission like CoRoT 
and Kepler. A current distribution of confirmed exoplanet^] can be seen in Figure 1.5 



(Schneider et al. 2011). 



There are various methods to detect an extrasolar planet. A detailed overview of all 



the methods can be found e. g. in Doyle (2008). 



1.2.1. Direct imaging 

The easiest way to detect an exoplanet would be the direct observation. This effective 
method has to deal with the problem that the light of the star is several magnitudes 
brighter than the planet. The brightness relation depends on two facts: 



The albedo of the planet. Since planets are not emitting light by themselves, except 
thermal emissions (Charbonneau et al. 2005), the only light source is the host star. A 
high albedo is required to measure reflected light from the planet. We may separate 
two classes of planets: rocky planets and Jupiter like planets. Taking Mercury as a 



template for a rocky planet one might expect a geometric albedo of 0.142 (Mallama 



Data taken from 



http://exoplanet.eu as of July 2011. 
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et al. 2002) and 0.52 for a Jupiter like planet]^]. The albedo of an extrasolar planet can 
be directly measured with the observation of the secondary transit a. k. a. occultation 



of the planet by its host star (Alonso et al. , 2009). We must be aware that our solar 



system is a rather bad template, and we might observe higher (e. g. Demory et al. , 2011 ) 



or lower albedos (e.g. Rowe et al. , 2008) 



The semi- major axis of the planet. The albedo of the planet has to be multiplied 

with the semi-major axis a, because the measured intensity of the planet varies with 
its distance from the host star. 

The direct observation can be improved by using a coronal mask, where the light of 
the host star is reduced, or using a nulling interferometer, where the light of the host 
star is cancelled out or choosing an appropriate wavelength, like the far infrared, where 
the planet is relatively bright. This has been successfully accomplished with the Hubble 



Space Telescope on Fomalhaut by Kalas et al. (2009), see Figure 1.7 



The biggest problem is the necessity for high resolution, and therefore for large aper- 
tures of interferometers. The next generation of observatories like the E-ELT or in- 
terferometric missions for the discovery of extrasolar planets are on the right way to 
overcome those limitations. 



1.2.2. The radial-velocity method 



Since the discovery of 51 Pegasi by Mayor et al. (1995), the discovery by radial- velocity 
(RV) has brought the majority of extrasolar planets. The mass of the extrasolar planet 
attracts the host star throughout the circular orbit and implies a sinusodial signal on 
the spectral lines of the host star. These lines are either blue- or red-shifted. 

The velocity signal V of the host star is made up by the host star's relative motion 
V in respect to the observer and the movement of its barycentre v*. It is given by 



V(t) = V, z + 



2TcaM p sin i 



(cos(#(£) + wqp) + ecoscuop) 



(M p + m*)pVT^ 

where a is the semi-major axis, Mp and M* are the masses of the planet and the star, % 



http : / / nssdc . gsf c . nasa . gov/ planetary/fact sheet / j upiterf act . html 
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is the inclination of the orbital plane, P is the orbital period, e is the orbital eccentricity, 
9 is the true anomaly and ojop is the angle between the pericenter and the orbital plane. 
The z-axis represents the vector from the observer to the star, wheras the x- and ?/-axes 
span the sky-plane. The expected amplitude 



.4 



RV 



2rraM p sin i 



(M p + M*)P^/T^ 



;i.s) 



can be calculated from Equation (1.4) 



Up to date all discoveries are verified by radial-velocity. Especially in combination 



with the transit method (see Section 1.2.7), the mass of the planet can be fixed. Having 
no other constraints on the planet, the radial- velocity suffers from the uncertainty of 
the inclination i implying an additional factor of sini on the mass of the planet. This 
means an amplification of 1.58 on average. 



Stellar activity puts a limit on the detection of exoplanets (Boisse et al. , 2011). Low 



mass planets can not be detected due to stellar RV jitter. For example the verification 
of Corot-7b took more than 80 nights on HARPS. The weak mass-determination of 



CoRoT-7b is caused by the very active host star (Hatzes et al. , 2011). 



1.2.3. Spectroscopic measurement of atmospheric features 

Another spectroscopic method similar to direct imaging would be the measurement of 
atmospheric species or molecules on the surface of the planet. Again, large apertures 
are required, to observe directly e.g. water or even chlorophyll on an illuminated planet. 

This topic is also related to transit spectroscopy, where the absorption of hydro- 
gen, sodium or water molecules can be detected. This kind of observation, additionally 



provides information on the evolution of the planet and especially its atmosphere (Lam 



mer et al. 2011). 



1.2.4. Astrometric method 

Since the star and the planet orbit a common center of mass, a perturbation in the stars 
motion can be detected with astrometric methods. The planet must be sufficiently large 
and the inclination close to zero to detect a variation in the stellar motion. This obeys 
the limit of detecting Jupiter class planets in very tight orbits. The astrometric method 
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has been probed by Lazorenko et al. (2011) on VB 10 as seen in Figure l.£ r 



The astrometric wobble w, where 

MpQp 



w 



MA 



1.6) 



and Mp is the mass of the planet, M* is the mass of the star observed at a distance d. 



1.2.5. Pulsar-timing 



Pulsars are known to emit periodic radio-pulses. This effect is caused by a strong 
magnetic field which poles are not in alignment with the rotation axis of the stellar 
remnant. Electromagnetic radiation is emitted, and the beam emerges into space. In 
rare occasions the pulsar is aligned in a way, that the beam hits the Earth. 

A massive body orbiting a pulsar causes variations of these radio emissions. Since 
habitable planets are not assumed around pulsars, they are no longer included in the 
exoplanet search. The first detection of this kind was a Jupiter-like planet orbiting the 



binary system PSR B 1620-26 detected by Backer et al. (1993). 



1.2.6. Detection due to Microlensing 

As a planetary system passes the line of sight to a distant star, the passing objects act 
as a gravitational lens according to their masses. The effect of the lens is, that the 
light is focused and hence intensified to the observer. These events occur only once per 
system but the duration of microlensing events is rather long. Due to the nature of the 
event, only Jupiter mass planets at several AU can be detected, which are definitely 
out of range for a confirmation by radial- velocity. 



1.2.7. Detection of Transiting Exoplanets 

We assume a planet that is orbiting a star. The inclination of the orbit is in the line 
of sight, so that the projected disk of the planet covers the star. This transit can be 
detected as a change in the light curve. The transits of Mercury or Venus, that can be 
observed from Earth in the solar system demonstrate this principle. To this day we are 
far from observing transits from Earth-sized planets. 

7 http : //www. nasa. gov/topics/universe/f eatures/C-VBlOb- 20090528 . html 
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Observing a solar eclipse is a special case of a transit, that is caused by the finite 
distance of the observer to the transiting object. As for stars we can always assume 
an infinite distance, that has no impact on the projected size of the planet on the host 
star's disk. We observe D ^> a, where the distance D is in the order of light-years and 
the semi-major axis a does not exceed several hundred AUs. A general exception from 
this rule are free floating planets, which are not discussed in the frame of this work. 



1.3. The Physics of Transits 



A detailed discussion on the transit parameters can be found e.g. in Seager and Mallen 



Ornelas 


(2003 


) or 


Winn 


(2010 



The depth of a transit AF is the relative change in flux 

' Rr 



AF 



R* 



;i.7) 



depending on the radius of the planet R p and its host star i?* as shown in Figure 



1.10 This is the simple geometric relation assuming an opaque planetary disk and 
uniform illuminating stellar disk. The planet itself does not emit any radiation either. 
Even in the far infrared the light emitted from a planet as thermal emission is several 
magnitudes fainter than the host star. 

The impact parameter b is defined as 



R* 



cost 



\ 



1 - v / AF) 2 - (Wt T ) 2 (l + 
1 - {t F /t T ) 2 



depends on the semi-major axis a and the radius of the star R* and the inclination i. 
In other words, the impact parameter is the distance between the transit path across 



the stellar disk and the center of the star as shown in Figure 1.10 It can be calculated 
from the measured lightcurve using the change in flux AF, the total transit duration tr 
and the duration of the flat part of the transit tp. High impact factors or inclinations 
respectively lead to a grazing transit. This type of transit is V-shaped and has to be 



distinguished from other blend scenarios as shown in Figure 1.4 
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The transit duration tx with an impact factor b = is 

t T = P S m- 1 ( R ^). (1.9) 



7r \ a 

By introducing 



I = ^ (R^ + R p ) 2 - a 2 cos 2 i, (1.10) 
we get 

Pr = P sin" 1 f - ] . (1.11) 
7T \a J 

The transit duration can be approximated under the assumptions R p <C R* <C a with 

tT~2R,^, (1.12) 

where G = 6.67300 x 10~ 11 m 3 kg~ 1 s~ 2 is the gravitational constant. 

The semi-major axis a depends on the measured period P using Kepler's third law 
a3= P>G(M. + M p)< (113) 

and can be simplified under the assumption M p -C M*. Other parameters can be 
directly derived from the lightcurve 

- = 2 -^4^L, (1.14) 
32P AF 3 / 4 

p* = — (1.15) 



A more detailed derivation of these parameters and can be found in Perryman (2011) 
and references therein. 

Additionally by the stellar mass-radius relation 

R* = kM*, (1.16) 



20 



1 . 3. The Physics of Transits 



and by taking the values k = 1 and x ~ 0.8 from Cox (2000) for main sequence stars 
we obtain 

v^f. (1.17) 

Po/ 

A less physical parameter is the probability of detecting a transit 

P = 0.0045^ x x 1 + eCQ t /2) - (1.18) 

It is assumed that the planets' orbits are distributed normally. Current calculations 
predict one planet for each star on the average or at least one planet among four solar 



like stars Borde et al. (2003) 



To simulate a transit numerically the following parameters have to be taken into 
account: 

• the radius of the star 

• the radius of the planet 

• the distance between the star and the planet 

• the distance between the observer and the star 

• the center to limb variation of the star's atmosphere 

• inclination of the planet's orbit 

Since not all parameters can be observed or measured respectively, we have to reduce 
them to something useful: 

• fractional radius of the planet 

• the impact factor of the planet 

• the center to limb variation of the star's atmosphere 

There are several levels of accuracy to simulate a transit of an exoplanet. The intensity 
I of the star can be normalized to 1.0. 
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1.3.1. Simple Transit Modelling 

The most simple approach separates the transit of a planet with radius R p across a star 
with radius -R* into two distinct cases: 

2nRl for d > R* 

(1.19) 

2nRl - 2%Rl else 

The variable d represents the distance between the center of the two discs. Needless to 
say this intensity can be normalized as stated above to: 

1 for d > i?* 

2nRi else 

A discontinuity arises at the point where d = R*. 



1.3.2. The Continuous Transit Model 

To provide a continuous function that approximates the observation, a geometric rep- 
resentation for this problem has to be found. 

We have to introduce some variables to calculate the overlapping disc segments A\ 
and A 2 : 



I = R p sin arccos 



R 2 — d 2 — i?p 
2dR n 



hi = R p - JRl - P, 



h 2 = R*~ yjR* - P, 

A 1 = f(R p ,h 1 ), 
A 2 = f(R*,h 2 ). 



;i.2i) 

;i.22) 
;i.23) 

1.24) 

;i.25) 



Where I ist the half length of the intersecting line and hi, h 2 are the heights of the 



stellar and the planetary segment, see Figure (1.11) for details. The formula for the 
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area of a disc segment with radius R p and height h is 



A(h, R p ) = i?p arcsin 



2R T 



(R p -h)J2hR p -h 2 . 



;i.26) 



We can now calculate the intensity 



I(d, R*, R p 



2ttRI 

2ttR 2 - 2nR 2 p 
irRl -A 1 -A 2 



. . . d > -R* + -Rp 
. . . d < .R* — .Rp 



7ri? 2 - ttRI - A 2 + A, ...d 2 >Rl + R 2 v 



;i.27) 



The function (1.27) is analytic in d within the whole interval [— oo, oo] for given R* and 



-R p . 



1.3.3. Stellar Limb Darkening 

To introduce the impact factor h where 0<6<-R*orO<&<l respectivly the distance 
can be expressed as 



d = Vx 2 + b 2 



(1.28) 



where x is a function of time t and can be assumed linear. Inclination exceeding 
this boundary will lead to partial eclipses. To approximate the observations we need to 



introduce stellar parameters, at least the center-to-limb variation ( Claret , 2000 ) which 
can be described by 



1(0,6) 2/3 

, : = - 1 + - cos 9 
7(0,0) 5 V 2 



;i.29) 



which can be derived from the Eddington approximation. Since we need a function 
not dependent of an angle but as a function of the radius or the distance, we modify 
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equation |1.29| to 



2 / 3 
1(d) = - I 1 + - cosarcsin(rf) 



1.30) 



Usually quadratic limb darkening is used where 



\x = cos I 



y/l - r 2 , < r < 1 



[1.31) 



The intensity / is then defined as 



/(r) = l- 7l (l- /U )- 72 (l- / i) s 



1.32) 



where 7i + 72 < 1 are the limb darkening coefficients. 



1.3.4. Eclipse of a Uniform Source 



An analytic description of planetary transits is given by (Mandel and Agol, 2002) 



where I want to summarize some considerations. The model describes a transit of a 
opaque spherical planet on a spherical star. We are interested in the relation between 
unobscured to obscured flux F e (p, z) = 1 — X e (p, z) with neglecting the limb-darkening 
effect, we get 



\ e (p,z) 



P 2 K + Ki 



4z 2 -(l+z 2 -p 2 ) 2 
4 



1 + p < Z 

|1 — p\ < Z <1 + p 
Z < 1 — p 
z < p — 1 



1.33) 



where d is the center-to-center distance, between the star with radius i?* and a planet 
with radius R p . The normalized separation of centers z = d/R* and p = R p /R* the 
size ratio with 



Ki = cos 



^1-p 2 + z 2 

2z 



1.34) 
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and 



k 



cos 



_! p 2 + Z 2 - 1 

2pz 



;i.35) 



1.4. Overview 



The different detection methods will be explained briefly in Section 1.2 We will discuss 
the transit detection method in detail and focus on the parameters that can be derived 
from transit observations of exoplanets. 

We will briefly describe the data reduction pipeline used for observations at the obser- 
vatory Lustbuhel. An overview of the most important and successful transit detection 



algorithms is given in Section 2.2 Different transit modelling codes including codes 



that initially were used to describe stellar binaries will be discussed in Section 2.4 



Finally the effects introduced by ground- and space-based observations are discussed 



in Section |2.7| Several active stars hosting transiting planets are discussed in Sec- 
tion 
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Figure 1.5.: Semi-major axis in Astronomical Units (1AU = 1.49598 x 10 n mj versus 
planetary mass in Jupiter masses (Mj np = 1.89 86 x 10 27 kg. The hori- 
zontal dashed-dotted lines represent the transit depth for an Earth- and 
Jupiter-sized planet, the dashed lines show the detection limits for ra- 
dial velocity at lm/s and lOm/s respectively. The vertical lines mark 
the maximum observable semi-major axis CoRoT (150 days) and Kepler 
(12 years lifetime assumed). The empty domain between a = 0.01 and 
a = 0.04 is called the Neptune void. 
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Figure 1.6.: The Doppler and astrometric measurements are complementary, depend- 



ing on the inclination angle as seen by the observer (adapted from Brown 



ffijOffy ). A transit can only be observed in a small cone with an inclination 



90. 
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Fomalhaut System 

Hubble Space Telescope • ACS/HRC 



NASA, ESA, and P. Kalas (University of California, Berkeley) STScl-PRC08-39a 



Figure 1.7.: Image of Formalhaut b. Credits: \NASA\ \ESA\ P. Kalas, J. Graham, 
E. Chiang, E. Kite (University of California, Berkeley), M. Clampin 
(NASA Goddard Space Flight Center), M. Fitzgerald (Lawrence Liver- 
more National Laboratory), and K. Stapelfeldt and J. Krist (NASA Jet 
Propulsion Laboratory) 
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<- 1999 Sept 16 
<- 2003 Juy 7 

£ : 2008 Sept ! 

5 • 



Figure 1.8.: The planet VBlOb induces a wobble of 6 mas on its central star. 




Figure 1.9.: A planet transiting its host star, with the (pri mary) transit and the sec- 
ondary transit (occulation) . Figure taken from Winn (2010). 
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Figure 1.10.: Sketch of a transiting planet shown in gray in front of a star in yellow. 

R* is the radius of the star and b is the impact factor. The total transit 
time tx is the time between t\ and £4. The "Hat part" tp is the time 
between t 2 and £3- 




Figure 1.11.: Geometrical illustration of a transiting object. The star with radius R* is 
partially covered by the planet R p . The intersecting circles are separated 
by the line I. The heights of each segment is hi and h 2 respectively. The 
size of the planet has been scaled for better visualization. 
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Transit detection in the frame of active stars does not only depend on fitting the right 
model to the data, it affects the whole process from data acquisition to the resulting 
model. 



The process of transit observation is described in Section |2.1| Since transit depths 
are in the domain of a few mmags, care has to be taken for observational planning as 
well as for the data-reduction. 

Usually one in thousand stars or even less have a transiting exoplanet. This implies 
a large amount of stars to be observed simultaneously to ensure a detection. This 
principle is carried out on space-based and on ground-based observatories. Transit 
detection algorithms that can handle the data gaps caused e. g. by the night-and-day 
cycle are necessary to sieve the large number of targets. A selection is presented in 



Section 2.2 including the de-facto standard by Kovacs et al. (2002). I present my own 
considerations on Section 12.61 

Various filtering methods are discussed to increase the signal-to-noise ratio in Sec- 



tion 3.2 Special care is taken on the analysis of periodic signals in Section 2.5 



2.1. Observation of Transiting Exoplanets at the 
Observatory Lustbuhel 

I want to give a short description, how transiting exoplanets are observed at the Ob- 
servatory Lustbuhel Graz, and describe in detail the data reduction pipeline. Both 
telescopes, the 40 cm Schmidt-Cassegrain and the 30 cm Zeiss refractor are capable of 
observing transits on brighter targets. 
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2.1.1. Target Selection 

The targets have to fulfil a set of criteria to be observed. Usually the target is selected 
from a list provided by the Exoplanet Transit Database (ETD r] 




• The transit depth must be at least 0.01 mag to ensure a detection of the transit 
at a reasonable signal-to-noise ratio. 

• The V magnitude must be brighter than 15 magnitudes. 

• The airmass of the object must be less than 1.8. 

• The lunar phase must be less than 50% and the angular distance between the 
target and the moon should be at least 90°. 

• The target should be in an azimuth range between 60 and 120 degrees to avoid 
regions with light-pollution near Graz. 

2.1.2. Observation 

The observation starts at least 15 minutes before the transit to measure the level of 
stellar activity out-of-transit. 

A set of at least 20 dark frames is taken after the observation with the same exposure 
time as the light frames. Also a corresponding set of bias frames is taken at the same 
temperature to ensure the reduction. Flat fields are usually taken from a library. 

2.1.3. Image Registration Techniques 

Inaccuracies in tracking the target as well as guiding errors lead to image displacements. 
Higher airmasses usually imply higher guiding errors. To measure the flux of a target 
star and compare it to comparison stars, a constant position of the target is mandatory 
throughout the whole image stack. This can be basically accomplished by two different 
methods, either by keeping the image coordinates in pixel constant, or by calibrating the 
image to calculate the absolute positions on the sky in right ascension and declination. 

J http : / / var2 . astro . cz/ETD/ index . php 
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Fourier Image Registration 

If we assume a displacement in x- and y-axes between an image i and a reference r, the 
easiest way would be to probe all possible displacements starting from (0, 0) to (w, h) 
where w is the width of the image and h is the height of the image. Actually the image 
shift can be assumed in the range [—w/2, w/2] and [—h/2, h/2]. To calculate the match 
between the two images % and r, we calculate \% — r\, which we want to minimize. Using 
the Fourier cross-correlation 

csrW-j'M), (2.i) 

between the both images i and r, where J 7 is the Fourier transform and J 7-1 its inverse 
and J 7 * the complex conjugate. The Fourier transform of a 2D function is calculated 
by 



T(u, v) = J J f(x, y)e- 2 ^ ux+ ^dxdy, (2.2) 
where u and v are spatial frequencies and its inverse 

T~ x (u, v) = f(x,y) = J J T{u, v)e 2i7T(ux+vy) dudv. (2.3) 
The shift theorem in Fourier-space 

X[f(x - a)} = e - 2iltua 7 (2.4) 



The coordinates of the maximum in the image-array c in Equation 2.1 marks the dis- 
placement. This registration is only accurate by integer pixel values, which anyway is 
flux conserving. If we want to register the images in sub-pixel accuracy, the images 
have to be scaled up. So an upscale by the factor 2 will result in 1/2-pixel accuracy. 
If displacement and rotation have to be corrected the Fourier-Mellin transform has 



to be applied. An appropriate algorithm is described by Reddy and Chatterji (1996). 



Astrometric Registration 



Stellar-like sources can be easily identified and extracted with S Extractor by Bertin 



and Arnouts (1996). Each star is fitted with a two-dimensional Gaussian, where its 
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photo-center can be determined numerically. This method provides sub-pixel accuracy 
down to 10 milliarcseconds. The limitations are given by the detector noise and by the 
distortion caused by imperfect optics. 

By matching the stars to a precise catalogue, the exact transformation from pixel 
coordinates to celestial coordinates can be determined. In this way the target star is 
located by right ascension and declination and the stellar flux can be either measured 
from the gaussian fit done by SExtractor or by applying aperture photometry using 



the DAOPhot routines (Landsman, 1995). In latter case the stellar flux is measured 
within a circular area and the average sky-flux that has to be subtracted is measured 
in an annulus. 



2.1.4. Data Reduction 

In the whole dataset the target star V and two comparison stars C and K are selected, 
both being known to be non-variable. One of the comparison stars should be dimmer, 
the other brighter than the target star. All three stars must be present throughout 
the whole dataset, otherwise the frames where one star is missing is omitted. The 
comparison stars should have a photometric precision better than 0.01 magnitudes. 
Usually the photometric magnitude is mapped to a R-, J- or K-magnitude in a catalogue 
regardless of the fact that the measurements were taken in the visual or R-band. 

For each frame the observed values of the comparison stars are mapped to the known 
magnitudes resulting in a correcting linear function, where the precise photometric 
value for the target star is calculated by interpolation. This could be in principal be 
extended to a whole set of comparison stars by applying a linear regression. 

For high precision, the comparison stars should be of the same color temperature as 
the target star. This is generally not the case due to a limited field-of-view. Therefore 
as a last step the data has to be corrected for influences caused by the airmass. 



2.2. Detection Algorithms 

Various methods for the detection of transit signals have been developed. Most of them 
require a more or less "clean" dataset or lightcurve. Usually the data contains gaps and 



other irregularities like jumps, hot pixels (see Section 2.7.1). The detection of transit 
like features can be achieved in different approaches: 
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Brute force The simplest idea would be, to scan the lightcurve for recurring transit- 
like shapes. The parameter space is spanned by periods, depths and different impact- 
factors. 



Fold and Search The lightcurve can be folded with different periods, and the transit 



signal can be found in the folded data (see Section 2.5.1). 



Correlate with Model Depending on the spectral type of the star, a model can be 
generated in respect of the limb-darkening. This model lightcurve can be correlated 
with the dataset. 

A description of a Monte Carlo Markov Chain (MCMC) analysis of the parameters 



derived from the lightcurve can be found in e. g. ( Burke et al. , 2007 ) . A good comparison 



of performance of the algorithms described in the following sections can be found in 



e.g. (Tingley 2003 Moutou et al. 2005). 



2.2.1. The BAST algorithm for transit detection 



The BAST (Berlin Automatic Search for Transits) algorithm (Renner et al. 2008) 
was developed in advance to the CoRoT mission. It could not be tested on real data 
but simulated lightcurves for CoRoT. The basic principle is the search for box shaped 
signals in pre-filtered lightcurves. The filtering includes normalization, application of a 
low-pass filter and correction for variability. 

After filtering, the algorithm performs a search for periodic signals. The performance 



is similar to BLS (Section 2.2.2). Since BAST was only tested on simulated data, the 



results on real data are degraded. Renner et al. (2008) claims to be able to find transit 



depth down to 0.01% corresponding to an Earth-sized planet around a solar-type star. 



2.2.2. A box-fitting algorithm in the search for periodic transits 
(BLS) 



The algorithm by Kovacs et al. (2002) performs a box-fitting to the lightcurve. The 



transit is here assumed as a box-shaped function with two levels H and L. The box- 
function is fitted to a folded time series at a trial period. The period is then permuted. 
The algorithm weights different levels L, giving higher levels a higher priority. 
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A signal-to-noise-ratio of 6 is required to ensure a significant transit detection. The 
BLS algorithm is the current standard algorithm performed for transit searches. Its 
performance can be improved by filtering the lightcurve prior to processing. BLS is 
computational costly on the downside. 



2.2.3. TRUFAS, a wavelet-based algorithm for the rapid detection 
of planetary transits 



The Trufas algorithm by Regulo et al. (2007) is based on a continuous wavelet trans- 



form of the lightcurve, which has been detrended in a first step. In a second step 
the period is searched within the wavelet transformation. At least three transits are 
required to ensure a proper detection. 

This algorithm had been developed on simulated CoRoT-data as the BAST al- 



gorithm. With Trufas the preprocessing is not crucial. Regulo et al. (2007) claims 
that the algorithm is well suited for large data volumes. 



2.3. Filtering Methods 



An attempt of filtering of lightcurves in the presence of stellar activity has been demon- 



strated on one of the CoRoTBlind Tests by Alapini and Aigrain (2009). Filtering is 



necessary to remove spurious signals from the lightcurve and therefore to increase the 
signal-to-noise ratio (SNR). 

In our case the signal is the transit signal and the noise is the standard deviation of 
the out-of-transit data. All measurements show intrinsically Gaussian noise as a first 
approximation. For faint stars the photon noise becomes relevant, which is proportional 



to yN } where N is the number of photons (see Weingrill et al 



2011 



Aigrain et al. 



2009). Usually lightcurves are poised with red noise caused by the finite observational 



length and transient stellar activity features. 

Usually the model-fit of the transit has to be better than 3a in the unfiltered folded 
data. In rare occasions a marginal detection of la is accepted, when the lightcurve is 



noisy or one or more transits are missing. Assuming a model \i for the transit (2.4) 
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with n x data-points in transit and performing a reduced chi-squared fit 




of n p parameters with a variation of out-of-transit. 

2.3.1. Median Filter 

The median filter is a low pass filter, where a datapoint Xi is replaced by the median of 
its surrounding data Xi-k, ■ ■ ■ , Xi + k- The length of the filter corresponds to the charac- 
teristic timescale r. Usually high frequency content with r < 0.5 days is removed with 
this filter. 

By subtracting the median filtered data from the unfiltered original, a high pass filter 
is realised. It works ideally to remove long term variations caused by stellar activity. 



2.3.2. Wiener Filtering 

The Wiener filter operates in Fourier space under the assumption that the noise N(f) 
can be estimated in a power spectrum. So the optimal Wiener filter W(f) for a measured 
signal S(f) is defined as 

W(f) = \ S iDl (2.6) 

A priori we have to estimate the noise background independent of the frequency /, so 
the noise is a constant background in the power spectrum. 

In the case of planetary transit searches the instrumental effects and the stellar activ- 
ity show generally higher levels at the lower frequency end of the power spectrum which 
is called red noise. Wiener filtering is not suitable for eliminating stellar activity since 
it fails to estimate the noise background. 



2.3.3. Butterworth Filter 

The Butterworth filter consists of a flat part that leaves the wanted frequencies un- 
changed and a slope of variable steepness to block the unwanted frequencies. The gain 
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G{uj) is defined as 

G 2 H = f\ 2n , (2.7) 

where n is the order of the filter determining the steepness of the slope, u c is the cutoff 
frequency sets in and Go is the gain at zero frequency, which is usually unity. 

This filter can be applied as a substitute to the median filter, as a high- and a low-pass 
filter. 



2.3.4. Savitzky-Golay Filter 



The filter by Savitzky and Golay (1964) is generally speaking a smoothing filter for 
spikes. For that purpose the adjacent points . . . , Xi-x, Xi + i, . . . , Xi + k of each data- 
point Xi is interpolated by a low order polynomial. The value of Xi is replaced by the 
interpolated polynomial value. 

As a digital finite impulse response filter, the data is usually convolved with the 
coefficients of the filter kernel. The filter works fast and removes single outliers in 
the data without affecting the shape of the transit. The filter is mostly used to re- 
move instrumental noise, caused e. g. by spurious outliers from the Southern- Atlantic- 
Anomaly (SAA). 



2.3.5. Pre-whitening 

Stellar signals that are generated by rotation or pulsation can be approximated with 
a small set of sinusoidal functions. By determining the most prominent frequencies 
including their harmonics, those rotation signals can be eliminated. The filter works 
iteratively by performing the following steps: 

1. Calculate the Fourier power spectrum and determine the maximum amplitude 
with the respective frequency. 

2. Fit the function A ■ cos(wt + <ft) with the amplitude A the frequency u and the 
phase 0. 

3. Evaluate the stop condition (e. g. remaining total frequency power) and repeat 
step 1. 
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Usually fitting a cosine function (Section A.2)is more robust than a sine function. 
The fitting by minimizing the least squares can be either done using the AMOEBA or 



the Levenberg-Marquardt (LM) (Press et al. , 1992) algorithm. The same algorithm can 
be used to separate the noise from the signal, by iterating the removal of frequencies 
unless the amplitudes are in the order of la. 

Pre-whitening is usually performed to clean radial-velocity measurements from un- 



wanted frequencies (e.g. Queloz et al. 2009). It can also be applied selectively to 
remove known instrumental frequencies. 



2.3.6. Transit Filter 

The basic idea for a specialized filter is that we are looking for transits. We know that 
the transits cover a distinct area in parameter space, periods P = 0.5...T/3 days, 
where T is the total length of observation, transit durations ^ = 1 ... 20 hours, depths 
AF = 0.01... 0.1%. 



1.000 



0.995 



0.990 



0.985 



I I I 

i '■ L 


I I I I I I 


r 

+ + 


i i i i i i 




+ + 






" " " "+" T " " 

+++++ 







i 






38 



39 



40 

time [hours] 



41 



42 




100 200 300 400 
freq. [mHz] 



500 



Figure 2.1.: left: Transit of a Jupiter-sized object with a duration of 2 hours and a 
depth of 0.1%. The sampling was 512 s. right: Power spectrum of the 
transit-shape. 



Using a Fourier filter, it is possible to filter out just the signal, which we want to 



detect. Using the low-pass filter as seen in the right Figure 2.1, we are able to filter out 



transit shapes that look like Figure 2.1 (left). The transit itself was synthesized using 
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the Mandel and Agol (2002) algorithm with limbdarkening coefficients u\ 
u 2 = 0.75. 



0.25 and 



This filter must be adapted to different transit durations, while the depth and the 
limb- darkening is uncritical. Usually a single filter functions covers a large volume in 
parameter space. 

2.3.7. The SARS Algorithm 

A filtering method that considers mostly instrumental effects has been developed for 



CoRoT by Ofir et al. (2010). It is based on the assumption that all lightcurves are 



affected by instrumental noise on the CCD in the same way, hence showing slight cor- 
relations. The algorithm tries to minimize those influences by generating a photometric 
reference. 

The SARS algorithm is necessary to reach the domain of ppm precision. Using this 
kind of cleaning procedure faint transits in the regime of 0.01% are visible. Nevertheless 
systematic errors can be introduced and transit signals are often below a 3a confidence 
level. 



2.4. Transit Modelling 

Many parameters of a transiting system can be derived by analysing the observed 
lightcurve and comparing it to a calculated model. Most of the codes are descendants 
from binary star modelling tools. The most prominent codes that have been applied 



to exoplanets' transits are Phoebe by Prsa and Zwitter (2007) and JKTEBOP by 



Southworth et al. (2009). 



Stellar limb darkening coefficients are usually taken from Sing (2010), which are avail- 



able from Tefjfrom 3500 to 50000 K including different surface gravities and metallicities. 
Other codes have been developed to test transit search algorithms, like the Uni- 



versal Transit Modeller (UTM) by Deeg (2009). One has to distinguish transit 



modelling from transit search algorithms. Transit models are usually specialized on 
many parameters like e.g. limb darkening coefficients, which themselves are dependent 
on the temperature of the host star. These are not suitable to detect transits in a set 
of lightcurves due to the mathematical complexity 
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Codes designed for binaries have the major drawback that the secondary is assumed 
to be a stellar object, which extrasolar planets are definitely not. As a workaround the 
surface temperature of the secondary is assumed to be 3500 K to satisfy the fitting of 
limb darkening coefficients. Nevertheless the results from these codes are questionable 
since other parameters like e.g. surface brightness, gravity brightening may falsify the 
output. 

Analytic lightcurves as a base for searching transiting planets have been discussed by 



Mandel and Agol (2002), which is today the de-facto standard for transit models. 



2.4.1. Phoebe 



The Phoebe program by Prsa and Zwitter (2007) uses the code by Wilson and Devin 



ney (1971) for binary stars. It is capable of deriving model parameters for different 



types of binares using radial-velocity and photometric measurements. In the case of 
transiting planets the model for a detached binary is used and the surface temperature 
of the planet has to be assumed 3500 K. 



This code has been applied by Poddany (2008) successfully on Wasp-2b. The big 



advantage of Phoebe, using a scripted version to automatically find the best model 
parameters, can not be used due to the high mass ratio between primary and secondary 
object. Mass, radius and inclination as well as the ephemeris of the mid-transit can be 
estimated using Phoebe. 

The program provides a graphical user interface (GUI), which greatly simplifies the 
manipulation of many parameters for fitting the parameters. The major drawback is the 
fact that the initial parameters for the fit have to be very close (within a few percent) 
to converge. Usually the fitting procedure halts after some steps without convergence. 



2.4.2. JKTEBOP 



The JKTEBOP code has been used by Southworth et al. (2009) to model several transits 



of transiting exoplanets. Similar to Phoebe (see Section 2.4.1), this code incorporates 



the Levenberg-Marquardt optimisation algorithm (Press et al. , 1992 ) and advanced limb 
darkening parameters. 



JKTEBOP requires several input parameters (2.1 ) as a first estimate, but works more 
robust than Phoebe and in an automated way. 
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Table 2.1.: Some input parameters for JKTEBOP used to calculate WASP-4. 



parameter 


value 


Sum of the radii 


0.21 


Ratio of the radii 


0.15 


Orbital inclination (deg) 


88.5 


Mass ratio of system 


0.0013 


LD law type for star A 


quad 


LD law type for star B 


lin 



2.4.3. Mandel & Agol Algorithm 

The core of the M&A algorithm is used in the Transit Analysis Package (TAP) 



by Gazak et al. (2011). With the given parameters from Table 2.2 we get the result 



seen in Figure 2.2 



Table 2.2.: Imput parameters for a transit model. For a explanetion of the limb dark- 
ening coefficients see Section 



1.3.3 



description 




parameter 


value 


p = -R p /-R* 




P 


0.1 


impact factor 




b 


0.5 


limb-darkening 


coefficient 


7i 


0.25 


limb-darkening 


coefficient 


72 


0.75 



The algorithm is computational intensive and is not suitable for transit detection. It 
is intended to perform transit fitting to derive the transit parameters p, b, 71 and 72. 

By fitting a trapezoid using the Levenberg-Marquardt least squares algorithm im- 
plemented in IDL (Markwardt, 2009), we get the parameters p = 0.107, with a dura- 
tion of tf = 1.92 and the impact factor b = 0.92, which is highly overestimated. The 
X 2 = 3.73 x 10~ 4 is sufficient. The uniform source is generally better fitted with p — 0.1, 
t T = 1.93 and b = 0.585 with af = 4.3 x 10~ 5 about one magnitude better. 



2.5. Analysis of periodic signals 

One initial problem is finding the main period of a signal. In transit detection, at 
least three transit signatures are required to assume a transiting object. In rare cases, 
where a high signal-to- noise ratio is evident, two transits will suffice, bearing the risk 
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-2-10 1 2 

stellar radii 



Figure 2.2.: Transit model after Mandel and Agol (2002) with parameters taken from 
Table 2.2, The solid red line includes the limb-darkening coefficients 
7i)72? whereas the dashed line is calculated for a uniform source. The 
residuals are plotted with an arbitrary offset at the same scale. 



of a binary or even a blend scenario. Having time-limited observation windows like for 
CoRoT 150 days at maximum, sometimes mono-transits are detected. These detections 
are candidates for follow-up observations to eventually detect further transits. 

We define a signal as period, if a signature occurs at equidistant time intervals. Hence, 
if we fold the signal at the length of the given interval, and rebin the signal we get a 
minimal root-mean- square, which is defined as 

rms = L 2 (x — x). (2.8) 



Rebinning 

Each measurement of the flux of a star requires a certain time of integration on the 
CCD. For CoRoT this is typically 512s but will be changed to 32s if a target is of 
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interest. For the detailed analysis, we are interested in constant time bins for the data. 
Hence we need to rearrange the data by conserving the flux. This process is called 
drizzling and initially intended for images 

By rebinning the data to larger time intervals we get the standard deviation (a) in 
each bin as an error estimation. This value can be later on used for fitting the data. If 
we rebin 32s cadence data to 512s, we get 16 measurements in each new bin. Data gaps 
smaller than 512s still leave at least one measurement in the bin, but show a higher a. 
By looking for transit signals which last three hours, we might even consider to rebin 
the data to e.g. one hour bins, which condenses the data by a factor of seven. However, 
this is dangerous, since a single outlier in the large bin may cause an alteration of the 
data, in a way, such a transit search algorithm may not be able to find a transit or even 
detects a false positive. 

The major disadvantage of the drizzling algorithm is the fact that bins with a high 
local derivative, e.g. in the ingress or egress of a transit show abnormal high cx-values. 



Folding 

For n data-points di and their time-indices £, we calculate 

t'i = UmodP, (2.9) 
where P is the period and which fulfils 

*min<4<^ (2.10) 

by sorting the times ascending and applying the sort to the data di, we get a new 
sorted dataset d\. This new dataset is now folded at the period P. If the period is not 
an integer multiple of the original time-interval we have unequal spaced data. Using 



the algorithm from the previous section 2.5, we rebin the data. The implementation 
can be found in Section IA.91 



http : //www. stsci . edu/f tp/ science/hdf /combination/drizzle . html 
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2.5.1. Finding the period by phase folding 

The easiest method to find the main period of a signal is by iteratively trying possible 
periods, where we fold the data and look for the "best fit" in our case, where the 
maximum of all (7j for each bin i is a minimum. We give an example on the binary 
SRaOl El 0359. The domain, where we look for a periodic signal can be restricted 
between 0.5 days and 11 days for a total run length of 33 days. The lower boundary is 
given by physical limits. Planets or in our case a secondary star would be within the 
Roche-lobe of the host star or the primary respectively. The Roche limit d is 

/ \ 1/3 

d = rj2^\ , (2.11) 

depending on the stellar radius r*, the density p a st of the star and the density p p of the 
planet. 

The higher boundary is a third of the observation period, which means that at least 
three transits can be observed. This limit is given by logical constraints: a single transit 
won't reveal a period, whereas two transits can not exclude a binary system. Hence for 
looking for an Earth-analogue at 1 AU we need to observe three years in total. 

We start with a rough estimation of the period, which can be deduced from the 
second local maximum of the autocorrelation of the signal or the maximum in the 
power-spectrum of the signal. A proper algorithm is e.g. 



1. initialize a function-minimizer e.g. AMOEBA with the estimated period P 



2. fold the data with the given period P as described in Section 2.5 



3. rebin the folded data to d\ 



4. calculate the max(crj) over all bins as the weight for the function. 



We reiterate from point 2. until a local minimum is found. 

For the case of the binary SRa01_El_0359, AMOEBA found 7.01408 days. 

Algorithms that minimize a function like AMOEBA or Levenberg-Marquardt however 
need good initial conditions. In our case, a good guess of the true period. 
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Figure 2.3.: The maximum of the RMS of each bin of the binary SRa01_El_0359. A 
total of 500 periods between 0.5 days and 11.5 days were probed, giving 
the best period of 7.02 days. 

2.5.2. Fourier Autocorrelation 

The Fourier cross-correlation of two functions / and g is defined as 



c = JF-\F{f)-r{g)), 



(2.12) 



where T is the Fourier transform, J 7 1 its inverse and T* is the complex conjugate. 
The Fourier auto- correlation of a function / is defined as 



a = JF- l {jF{f).r{f)), 



(2.13) 



as the cross-correlation of the function with itself (see also Section A. 3). 



By performing a sliding windowed auto-correlation over the data with the length twice 
a characteristic timescale (e.g. rotation period or orbit period), transient features can 
be followed. This method is superior in using a sliding Fourier window, since it "locks" 
on the most prominent feature without any presumptions. 
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Sliding auto-correlation may be applied to detect timing variations. It operates as a 
wavelet-filter in a specific frequency. 



2.6. Description of a Transit Detection Algorithm 



We want to give a short description of a transits detection algorithm, that has been 



developed in the frame of the CoRoT(Baglin et al. 2006) detection pipeline. After an 



observation run of CoRoT the raw (NO) data is undergoing two pipelines (Nl, N2) where 
the invalid data-points are flagged, e.g. for hot pixels and for crossings of the satellite 
over the South Atlantic Anomaly (SAA). The SAA causes a semi-periodic alteration of 
the on-board electronics and therefore impacts the data. 

The CoRoT data can be downloaded from the CoRoT Data Center at IAS archived 
Each run in separated between chromatic and monochromatic targets. The chromatic 



targets are usually the bright targets (Figure 2.6). 
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Figure 2.4.: Distribution of chromatic targets over magnitude for the CoRoT run 
LRaOl . The gray bars show the amount of monochromatic targets. 

In front of the the CCDs El and E2 for the exo-fields a grism has been inserted to 
spread the point-spread-function (PSF) in its spectral colors. This enlarges the PSF 
to several arcseconds. The image of the PSF on the CCD is called imagette, as seen in 

3 http : / / idoc-corot . ias . u-psud.fr/ 
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Figure 2.6 By selecting distinct columns in the imagette, the color information can be 



extracted. On fainter stars the whole intensity of an imagette is integrated. 




Figure 2.5.: Imagette of the CoRoT target #0102728404. The coordinates are given 
in absolute pixel of the CCD. The stars * mark the position of other 
CoRoT targets, the crosses x indicate faint background stars and the 
plus-signs represent the true position on the sky. Imagette taken from 



Exo-Dat (Deleuil et al, 2009). 



2.6.1. Data Preparation 

Usually a CoRoT-lightcurve has a cadence of 512 seconds. The on-board cadence is 
32 seconds which is kept for interesting targets and candidates in the so called Alarm- 
mode. Due to bandwidth limitations in the satellite downlink, only a few targets 
are available in high cadence. The onboard electronics downsamples the 32s-data to 
512s. The timing of each datapoint is given in CoRoT Julian Date (CoRoTJD) and 
heliocentric Julian Date (HelJD). To use an absolute Julian date, we have to add a 
constant 

JD = HelJD + 2451545.0 (2.14) 
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where the constant corresponds to 1 January 2000. The total number of datapoints n 
in a dataset is calculated by 

n = [168.75(maxt - mint)] , (2.15) 

where the constant is used to convert to samples per day. 

For detection purposes the signal-to-noise ratio (snr) is at highest for the integrated 
flux combining all three color channels. To minimize the noise all data is oversampled to 



512 seconds. This is done using the drizzling algorithm (see Section A.l and Section 2.5) 
which bins the data according to its timestamps. Then the mean of each bin is calculated 
and the rms in each bin is used as an error estimate. A better way to integrate one bin 
would be taking the median, which produces a better snr but the noise statistics differ 
from the on-board processing of CoRoT. 

Since the drizzling may produce empty bins resulting from the flagged data, this 
data has to be substituted for following algorithms, that assume equidistant data- 
points. For that purpose the missing data is simply linearly interpolated. The linear 
interpolation has the major advantage, that the power spectrum is not severely affected 
by the introduction of the artificial data-points, whereas other interpolation methods 
like quadratic, or spline interpolation induce several artificial frequencies. The optimal 
interpolation function would be the Fourier interpolation (for an implementation see 



Section A. 5) which minimizes the total energy in the power spectrum. This however is 



computational expensive. 

After interpolation the mean flux of the lightcurve is normalized to 1. This is to 
overcome numerical problems, since using the original flux values which are in the 
order of 10 5 . An alternate method is the conversion to logarithmic values. 

f' = -2.5 Xl \og 10 f + x 2 , (2.16) 

where x\ = 0.9674492 and 22 = 26.3424511. The value of 22 marks the so called zero 



point. Mazeh et al. (2009) has shown that this is however dependent on instrumental 



effects. Both values were determined by applying a linear fit to the magnitudes of all 



target-stars of the run LRaOl (Figure 2.6) 



With x\ ~ 1 and £2 in the order of the galactic background the values look reasonable. 
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Figure 2.6.: Mean R magnitude versus counts in a logarithmic scale for the run LRaOl. 

The linear correlation is evident. Chromatic targets are shown in red, 
where the flux from all three channels was summed up. 

2.6.2. Detrending 

Most of the light-curves show a clear more or less linear trend, that can be caused 
either by long term stellar activity or by instrumental effects. Especially for giants, 
the rotational periods are in the order of the length of one observation run. The most 
puzzling instrumental effect on the long runs of CoRoT is the aberration caused by the 
satellite's motion along the earth's orbit. The stellar positions on the CCD move in 
respect to the guide sensor, but as the imagette remains at a constant position, the flux 
varies between the color-channels. 

Generally a second order polynomial in the form 

f(t)=x + x 1 f(t) + x 2 f{t) 2 (2.17) 
can be fitted and subtracted from the light-curve, where f(t) is the flux. The fit 
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must be applied on original data-points and not on those added by interpolation (see 
Section 2.6.1 ). 

In rare cases, where the modulation by the stellar rotation is evident, a fit with a 
sinusoidal is applicable 



fit) 



a cos 



( 



2-Kt 



\ T 



+ +k, 



(2.18) 



with the period r, a phase shift 0, an amplitude a and an offset k. This fit can also 



be iteratively used for pre-whitening (Section 2.3.5). Using the cosine instead of a 
sine function has numerical advantages to avoid sinO = 0, where the amplitude is 
undetermined. 



2.6.3. Filtering 

There are various approaches to separate the transit signal from any other signal caused 
by stellar activity or instrumental influences. For a run of 150 days we expect the transit 
signal to be shorter than 12 hours for a dwarf star. Therefore we apply a median filter 
that operates as a high-pass filter on the data. The width of the filter is 84 data-points 
where the median for a sliding window is calculated which contains 84 measurements. 
As a precondition, the measurements have to be equidistant in time. Otherwise, if large 
gaps occur, the median filter will be broader than necessary. 



To remove spikes and single outliers the filter by Savitzky and Golay (1964) is very 
efficient. Outliers will be replaced by a polynomial evaluated through the adjacent 
values. Usually 12 values (six before and six after the data-point) are interpolated by 
a second order polynomial. The filtering is applied by convolution of the filter kernel 
with the dataset. 

As a last optional filter stage, a low-pass Butterworth filter can be applied. Other 
filtering techniques are application of a Gaussian wavelet filter or a filter using the 



Haar-wavelet. An example of a wavelet representation can be seen in Figure 2.7 A 
filter would cut a horizontal slice to isolate the transits. 

Wavelets in principal are a powerful tool to remove noise. One should keep in mind 
that any filtering also deforms the transit signals, sometimes to an unknown extend. 
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Figure 2.7.: Wavelet representation of CoRoT-6. The power of the signal is repres- 
ented by differet colors. The wavelets in the domain from 0.05 to 0.2 
represent the transit signatures. 



2.6.4. Transit Detection 

After the preparation and the filtering process the main signal that remains is (hope- 
fully) transits with minor residual noise. The last task is the detection of transit-like 



signals. The most prominent parameter is the transit period (Section 1.3) that can be 



accomplished by calculating maximum of the autocorrelation function (Section 2.5). 



At this stage it is sufficient to fit a trapezoid or a box-function to estimate the transit 
parameters like depth, duration and the impact parameter. 



2.7. Separation of Stellar Activity from Instrumental 
Effects 

While instrumental effects on ground based observations are not an issue, they have 
a big impact on space based missions. Most observational and instrumental effects on 
ground are the duty cycle of observation. Photometric measurements are interrupted 
by daylight, which imprints a spurious frequency signal corresponding to one per day or 
11.5741 |j.Hz. Other effects are the changing airmass during the observation, different 



52 



2. 7. Separation of Stellar Activity from Instrumental Effects 



seeing conditions and sporadic cosmic ray hits. 

Space is a rougher environment for doing observations. The radiation level is much 
higher causing pixel alterations by cosmic ray hits as the main source of disturbances. 
Other variations derive from the satellites orbit and the altitude control. 



2.7.1. Instrumental Effects 

Instrumental effects can affect the whole image or specific pixels. We might separate 
effects on their impact on the measurements. For the whole discussion in this section 
we focus on space photometry, since ground based photometry suffers from other limit- 
ations like observational gaps due to weather conditions and or the day and night cycle. 



Even for observatories located at the South Pole (Strassmeier et al. 2008) there are 
limits due to the length of the Arctic night and the stability of the conditions as well 
as local seeing conditions. 
Most important effects are 

• Spacecraft jitter 

• thermal noise 

• readout noise 

• cosmic ray hits 

We must not forget optical effects either wanted as in CoRoT, where a prism extends 
the PSF to obtain color information, a Fresnel mask like in MOST to increase the 
signal-to-noise ratio by expanding the PSF. We also have unwanted or unavoidable 
effects like diffraction spikes from the spider holding the secondary mirror. The last 
one has a mayor impact on multi target photometry missions like Kepler and CoRoT. 

Bright stars in the observed field cause diffraction spikes affecting adjacent columns 
and rows of the star. This can lead to weird results where multiple targets share the 
identical period and epoch if the diluting star is a binary. Even the CCD-readout may 
produce correlated events along the readout direction. 



Jitter 

Jitter is caused by the spacecraft pointing. Since a satellite is a free floating object, 
its pointing has to be maintained by closed loop positioning system based on a star 
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finder and reaction wheels. The pointing accuracy is in the order below one arcsecond. 
However the centroid of a star is slightly moved on the CCD where pixels do not have 
the same sensitivity. 

Thermal noise 

The thermal noise is almost neglectable as it adds uniformly signal on the detector. 
Usually the detector is kept at a constant temperature. Ground based observations may 
suffer from temperature changes and introducing noise at different levels. The thermal 
noise is usually removed using a so-called dark-frame. It has the same exposure time 
as the image, but with the shutter closed. 

Dark frames also include cosmic ray hits. Depending on the level of accuracy the 
number of dark frames is usually in the order of 20 or even higher for the mmag regime. 
For space based observatories dark-frames are taken in the calibration & validation 
phase prior to the scientific part of the mission. 

Readout noise 

With todays electronics the readout noise is negligible. Generally the pixel-to-pixel vari- 
ations have a larger impact on the data. Occasionally we can speak of pixel correlated 
noise where the signal is correlated with the position of the target on the CCD. 

Cosmic ray hits 

Cosmic ray hits cause a local alteration in the CCD with different results. The most 
prominent effect is a sudden rise combined with a slow decay back to the previous level. 
The effect of the cosmic ray hit can be seen in Figure |2.8 

Unfortunately cosmic ray hits have the same shape as stellar flare, the only difference 
is its time scale. Usually the pixel alteration of a cosmic lasts longer than one day, 
whereas a flare lasts only hours. Also flares affect all color channels almost equally. 

Hot pixels 

Pixels on the detector with a constant (high) saturation level or a non-linear behaviour 
are called hot pixels. If they occur on the border of a photometric mask, they will show 
abnormal activity due to spacecraft jitter which is often referred as telegraphic noise. 
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Figure 2.8.: Cosmic ray hit in the blue channel of SRa02_El_0241 at CoRoTJD 
3225.62. The other two channels remain unharmed. The linear trend 
(dashed line) shows an offset by almost 0.1 mag after the hit, which 
indicates a permanent alteration of the CCD's physical structure. A 
typical exponential decay is visible. 

Telegraphic noise can occur quasi-periodic or even periodic and might be misinter- 
preted as a transit signal. At low signal-to-noise ratios, this is hard to distinguish. 
We can categorize the effects in groups: 

• affecting the whole image in the same manner by means of ADU and noise stat- 
istics. 

• signal depending on the pixel coordinates 

• random events like cosmic ray hits 

Image affecting noise 

The first one is also referred to as bias and can be corrected using calibration frames 
acquired by the spacecraft. But usually this frames are only taken at the calibration 
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Figure 2.9.: We see telegraphic noise occurring in the blue channel of the CoRoT- 
lightcurve SRa02_E1 0381 . The green channel remains unharmed. The 
jumps have a depth of 30 mmag. Spacecraft jitter causes constant jumps 
occurring quasi-periodic. In this case the photometric mask was moved 
over a hot pixel. 




phase prior to the main mission. Calibration frames are not recorded during the ongoing 
mission since observation time is valuable on a spacecraft. The background signal is 
subtracted in die CoRoT NO-pipeline. 



Pixel correlated noise 



The second can be quantified using autocorrelation between stellar targets or correlation 



between pixel locations (stellar coordinates) and noise properties (Section 2.7.2). The 
correlation may show itself in a small area or cluster of pixels on the CCD and therefore 
on stellar coordinates or correlate to the distance of a coordinate. 
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Figure 2.10.: We see telegraphic noise occurring in the blue channel of the CoRoT 
lightcurve SRa02_E1 1648. The green and red channel are not affected. 
The jumps have a depth of 60 mmag. This lightcurve can not be repaired 
for the search for transit signals. 

Random pixel affecting events 

Events like cosmic ray hits cause a sudden change in the pixel ADU value. Whereas 
hot pixels cause a (almost) constant value. 



2.7.2. How to quantify noise 

Usually noise is identified in the frequency spectrum as a constant background con- 
tributing to all frequencies equally. This is so called white noise. It is easily identified 



and removed by Wiener filtering (Figure 2.11) or a Butterworth filter operating as a 
low-pass filter. 

Other quantification methods are e.g. the mean absolute derivative (mad) or median 
absolute derivative or the rms which exceeds a specific value to locally identify noisy 
events. Generally it is not easy to clearly identify noise. 

Stellar activity reveals itself as red noise since activity is always limited to certain 
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Figure 2.11.: Power spectrum of CoRoT-6 with signal and noise levels estimated with 
a polynomial. The signal is indicated with a dashed line, the noise is 
indicated by a dotted line. Both signal plus noise are represented by the 
dashed dotted line. We see that the signal is underestimated, while the 
level of noise is somewhat accurate. 

timescales like rotation periods, spot-lifetimes, granulation lifetimes, .... This amplifies 
low frequencies in different levels. 

2.7.3. Analysis of Stellar Signals 

Stellar variability can be separated into two main classes: 

• intrinsic stellar variability, e.g. pulsation, flares, Novae 

• extrinsic stellar variability, e.g. binaries, transits of exoplanets, stellar rotation 
(spots). 

Pulsation and rotation signatures are easily identified using either a Fourier spectrum 
or the autocorrelation function. Usually the latter is more robust. A sliding autocorrel- 
ation function with a certain window size (e. g. twice the stellar rotational period) can 
be used to observe transient features like stellar spots. The lifetime of spots and their 
appearance and disappearance causes a change in phase that might result in cancelling 
out the frequency of the true stellar rotation. 
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A semi-automatic approach to identify stellar variability classes in the HR-diagram 



has been developed and applied by Debosscher et al. (2009) on the CoRoT data. 

Other non-periodic phenomena like flares could be identified by using a model which 
is cross-correlated with the data. 
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3.1. Instrumental Noise 



The CoRoT telescope was equipped with a unprecedented baffling system to avoid 



parasitic stray-light in the optical path (Auvergne et al. , 2009). Nevertheless the orbital 



cycle of CoRoT and its disturbance in the SAA is visible in the data, which can be seen 



in Figure 3.1 



Though CoRoT is an exceptional example, the orbital cycle is not visible in general. 
It has been avoided by Kepler-mission in an Earth-trailing orbit. Smaller photometric 
satellites like MOST are more affected. 



3.2. Filtering 

A variety of filters for enhancing the signal-to-noise ratio in lightcurves was discussed 
in Section . As stated earlier filtering has to be applied carefully without harming the 



transit-signal. A result of the filer-pipeline is shown in Figure 3.2 

A planetary transit in a five day orbit was simulated with the transit duration in 
hours. The original sampling was 32 seconds corresponding to CoRoT's high cadence 
sampling rate. The original series consisted of 13 500 points. Generally the filtering 
process does not depend on the shape of the transit itself, so an 'average' transit was 
chosen with p — 0.1, b — 0.5 and 71 = 0.75, 72 = 0. The data was downsampled to 
512 seconds leaving 844 points. The general shape of the transit remains unharmed. A 
slight shift appears since the datapoints are rebinned to the beginning of the sampling 
time. 

The median high-pass filter was applied with a width of 84 datapoints, which did 
not change the transit at all. Using the median filter as a low-pass, the flat part of the 



transit was a little bit more enhanced (dark green in Figure 3.2). 
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Figure 3.1.: High frequency part of the power spectrum of a F9V star CoRoT 
LRa01_E2 0205. The orbital frequency of 0.323 mHz and its harmonics 
are indicated by dashed vertical lines. 
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Figure 3.2.: A simulated transit (black solid) and the results of the filtering process. 

The red crosses indicate the 512 second oversampled data, green after a 
media n filter was applied and blue after applying \Savitzky and Golay 
{196$. 
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At a last stage the filter by Savitzky and Golay (1964) was applied. This filter is 
essentially to remove outliers that have not been flagged as invalid datapoints in the 
pipeline before. This filter has a severe impact on the transit transforming it into a 
Gaussian. Ingress and egress are slightly enhanced. 



3.3. Solar Activity 



The sun may act as a template for stellar activity, but we have to keep in mind that the 
sun, as we see and measure it now represents a star at the age of 4.5 Gyr and a spectral 
type G2V. It is difficult to scale solar activity on photometric levels to stellar activity. 



An attempt so simulate transits at solar like stars has been conducted by Carpano and 



Fridlund (2008) 



Strong flare events that radiate in the continuum are called "white-light flares" in 
contradiction to the usual flare event that can be observed in Ha. Even moderate 



C-class flares contribute to the TSI ( Kretzschmar 2011) 



The Haloween event on 28 October 2003 was an exceptional event with an X17-class 
flare. The flare, that was also visible in the continuum increased the TSI for 265 ppm 



(Kopp et al. , 2004 Woods and Kopp 2005) which can be seen in Figure 3.2 L 



3.3.1. Solar Irradiance 



The Total Solar Irradiance (TSI) varies with the solar magnetic cycle (Figure 3.4). 



The output variation correlates with the total output and the phase of the solar cycle 



(Figure 3.5). 



Space missions like CoRoT and Kepler don't have observational data that covers a 
full stellar cycle, other photometric missions like e. g. Hipparcos do not have sufficient 
datapoints to identify a periodic behaviour. Since the TSI is available with a cadence 
of 6 hours, the CoRoT-data can be resampled to be comparable with the solar activity. 
Since LRaOl is one of the longest runs observed by CoRoT, this run with its 150 days 
is the ideal testbed to compare to solar irradiance. 



1 Image source: http : / / spot . Colorado . edu/~koppg/TSI/ 
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Figure 3.3.: First white-light Hare observed with a TSI (red diamonds) instru- 
ment. The NOAA/GOES X-ray flux is shown fo r comparison (blue 



squares). (Kopp et al. 



2004 



Woods and Kopp, 2005) 
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Figure 3.4.: Total Solar Irradiance measured with the TIM/SORCE instrument with 
a 6 hour cadence. The data covers a little less than a solar cycle. High 
activity signals correspond to a higher output of the sun. The largest 
event in 2003 was caused by a solar spot 



66 



3.3. Solar Activity 




1360.6 1360.8 1361.0 1361.2 1361.4 
TSI [W/m~2] 



Figure 3.5.: Correlation between the Total Solar Irradiance measured with the 
TIM/SORCE instrument and its variance over a timescale of 150 days 
(blue) and over 30 days (red). The lengths have been chosen in accord- 
ance with the CoRoT observational run lengths. 



3.3.2. Solar Rotation 



The mean synodic rotation period of solar sunspots is 27.2753 days, which was defined 



by Carrington (Stix, 2002, p. 278). The rotational signal is not clearly evident in the 
Fourier spectrum of the TSI. This fact is caused by the lifetime and the distribution of 
the solar spots. 



3.3.3. Solar Oscillations 



Solar low-degree p-modes are centered around 3mHz. Their energy does not exceed 



2.5 ppm (Finsterle and Frohlich, 2001). Generally they can be observed in the domain 



from 2.1 mHz to 3.9 mHz with lowest power at the limits between 0.3 ppm and 0.1 ppm 
flStbcl [20021. 
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Figure 3.6.: Power spectrum of a F9V star CoRoT LRa01_E2_0205. The spectral 
power decreases from low frequencies to high frequencies. A small peak 
at the one day period is visible, as well as stronger peaks for the orbital 
period at 0.323 mHz and its harmonics (see also Ballot et al. , 2011). The 
granulation noise is visible between 0.1 and lmHz. The p-mode oscilla- 
tions are centred around 4 mHz. 



3.4. Stellar Activity 



3.4.1. Stellar noise 



A detailed discussion of noise properties in the CoRoT data is discussed by |Aigrain 



et al. (2009), where the variability for some observed targets reaches 0.5 mmag for 
hours timescales. This behaviour can be correlated to the activity of giant stars. On 
the other side dwarf stars show median correlated noise levels over 2 hours that reach 
0.05 mmag at the bright end. 

The noise level of the R-magnitude of dwarfs increases with the length of the mission, 
which is expected from sensor degradation. In general we observe red noise. An example 



of red noise is shown in Figure 3.6 
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In order to compare the stellar activity levels with the activity of the sun, the long 
run LRaOl was investigated. The rms levels of this run are comparable with other runs, 
with a slight increase caused by sensor degradation. 

The light-curves for the run were rebinned to 6 hours using the drizzle algorithm 
(Section sec:rebinning). The relative average standard deviation was calculated, that is 
naturally magnitude dependent since the photon noise contributes stronger on weaker 
fluxes. To get magnitude independent measurements the level of photon noise (~ 
1/J(N)) was subtracted. 

To get the noise properties for each spectral class the stars were grouped by their 
color temperature, putting five subclasses in each group. The minimum, median and 
maximum rms in each bin was calculated with the results shown in Table 13.11 The 
minimum noise level is 144 ± 38 ppm for all classes which contains granular noise and 
a not definable instrumental noise. This might explain the higher levels for late type 
stars, which are at the fainter end. 



Table 3.1.: Spectral types of dwarf stars taken from Cox (2000) with their effective 
temperatures. The minimum and median rms are given in ppm, the max- 
imum rms in percent. The values for the spectraltype K is poised by binary 
stars. 



Spectral Type T e g min. rms me 


dian rms max. rms 


A5V 


8180 


123 


268 


1.62 


FOV 


7300 


101 


394 


1.69 


F5V 


6 650 


92 


526 


1.95 


GOV 


5 940 


121 


465 


6.79 


G2V 


5 790 


145 


450 


1.46 


G5V 


5 560 


104 


349 


2.38 


KOV 


5150 


160 






K5V 


4410 


182 






MOV 


3 840 


192 


420 


2.32 


M2V 


3520 


193 


213 


3.37 


M5V 


3170 


175 


310 


1.28 



The mean activity level of all spectral types is 377 ± 100 ppm, which is more than 



the average sun at 244.13 ppm, which is indicated in Figure 3.7 
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Figure 3.7.: Stellar noise levels within 6 hours on the CoRoT run LRaOl. The solar 
noise level is shown in red indicating the minimum and maximum noise 
signal during the solar cycle. 
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Figure 3.8.: Periods in cycles-per-day (cpd) versus amplitude of all stars in the CoRoT 
run LRaOl. The red line corresponds to a transit depth of Jupiter, the 
green line corresponds to a Neptune and the blue line to an Earth-sized 
transit depth. 



3.4.2. Stellar Rotation 

The long run of LRaOl with its 150 days provides an observational window of sufficient 



length to observe stellar rotation signals. In Figure 3.8 all stars (including giants) were 
probed for ther dominant frequency in the Fourier spectrum which is assumed to be 
either from rotational or pulsational origin. 

We see a pileup of rotations/pulsations for more than 3 days (~ 0.3cpds). Generally 
we would expect the amount of pulsators among these stars to be rather low. 



3.4.3. Stellar Oscillations 



Stellar oscillations of a F9V star can be seen in Figure 3J3 We have to clarify that these 
measurements are form the so called exo-field of CoRoT which is not intended to be 
used for precise measurements of stellar oscillations. Nevertheless a sampling frequency 
of 31.25 niHz is sufficiently high to show first hints of stellar oscillation signatures. 
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Figure 3.9.: Partial power spectrum of a F9V star CoRoT LRa01_E2_0205. The 
spectral power decreases from low frequencies to high frequencies. A 
small peak at the one day period is visible, as well as stronger peaks 



for the orbital period at 0.323 mHz and its harmonics (see also Ballot 



et al. 2011). The granulation noise is visible between 0.1 and lmHz. 



The p-mode oscillations are centered around 4 mHz 
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3.5. Exoplanets around Active Stars 

In this section, case-studies on several discovered exoplanets that transit active stars 
are performed, as well as a hypothetical body with two Earth-radii transiting the sun 
will be investigated. 



3.5.1. A Super-Earth transiting the sun 

The sun is the only star so far that can be spatially resolved. It is hence a valuable 
target to perform several simulations to gain knowledge about surface features. High 
resolution images from the NASA SDO data archive^] are used as the image source. 
The continuum wavelength of 450 nm provides a substitute for a broadband white light 



image. The daily image of 6 June 2011 as shown in Figure 3.10 was chosen for the 
simulations. The solar image shows moderate activities. 



Figure 3.10.: Solar image form 6 June 2011 taken with the SDO/AIA instrument at a 
wavelength of 450 nm. Four different quadrants are visible in the image, 
which might be caused by instrumental effects. 

The diameter of the sun was measured from the image with 3145.45 pixels. The 
total image has a height and a width of 4096 x 4096 pixels. The radius in pixels of the 

http : / / sdowww . lmsal . com/ 
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Figure 3.11. 



Simulation of a transiting planet with two Earth radii. The black sold 
line is the simulated transit measurement. The red line is the fit of a 
model by Mandel and Agol (2002). The solid green line are the residuals 
of the fit. 



transiting Super-Earth with 2i?Earth was determined 28.81 pixels. All calculations and 
the transit-simulation were performed in subpixel accuracy. 

For the transit-simulation itself, the planet was shifted pixel-by-pixel over the solar 
image, with the integral flux being measured and normalized with the integral flux of 



the empty image as seen in Figure 3.10 The result of this simulation can be seen in 



Figure 3.11 



The fit was performed using the Levenberg-Marquard algorithm of the Mandel and 



Agol (2002) algorithm. The parameters of the fit and its results can be seen in Table 3.2 



Table 3.2.: Parameters for fitting the model by Mandel and Agol (2002). The initial 



values represent the a priori parameters for the simulation. The limb 



darkening parameters were taken from Stix (2002, p. 162). The errors are 
the reduced \ 2 values derived from the fit. 



parameter 


initial value 


fit 


error 


b 


0.0 


-0.00147472 


0.011713614 


P 


0.068 


0.0580131 


1.4103004e-005 


7i 


0.6 


0.890749 


0.0042045146 


72 


0.0 


-0.165174 


0.0066625388 
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Figure 3.12.: Enlarged residuals from Figure 3.11 The diameter of the planet is rep- 
resented of the vertical dashed lines. The center of the star is indicated 
by the vertical dashed dotted line. 



If we constrain the limb-darkening parameters to be within [0 ... 1] and also the 
impact parameter to be of positive value, we get b ~ ± 0.015, p = 0.0580248 ± 
0.00001666, 7i = 0.79096 ± 0.00147 and 72 = 0. The limb darkening is far from 
the theoretical value of 71 = 0.6 given in Stix (2002). The initial value for p was 



estimated from y5F which is slightly overestimating, due to the contribution of the 



limb darkening. A detailed view of the residuals is shown in Figure 3.12 



Interestingly we see out-of-transit variations, especially before the transit at i£* ~ 
— 1.1. This feature can be explained by instrumental effects, e.g. insufficient baffling, 
or we see traces of the K corona, which is also reasonable. The major error results from 
the imperfect limb darkening function that is unable to describe the limb correctly. 
The second largest error results from the quadratic limb darkening model which is 
also insufficient. A fifth degree polynomial is necessary to describe the limb darkening 
correctly. What remains are traces of active regions between 0.6 • • • — 0.3 and 

i?* ~ 0.05 . . . 0.5. The asymmetry is caused by instrumental effects, possibly a gradient 
in the CCD sensitivity or the filter bandwidth. A non-linear response of the CCD might 
also cause a wrong determination of the limb darkening coefficients. 



3.5.2. Stellar Activity on C0R0T-6 



The transiting extrasolar planet CoRoT-6b (Fridlund et al., 2010) is an interesting 



case. Especially the discrepancy between the photometrically active lightcurve and the 
missing activity signal in the spectroscopically measured Ca II H&K lines, proves the 
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Figure 3.13.: Lightcurve of CoRoT-6 with the transits of CoRoT-6b. In the last 
quarter the stellar activity signal is in the order of the transit depth. 



necessity to investigate transit lightcurves. The 'warm Jupiter' planet orbits the active 
star in a 8.9 day orbit in a distance of 0.08 AU. The host star is of solar type (F9V) 
and has an estimated age of 3.0 Gyr. 



The level of noise is rather high in the lightcurve of CoRoT-6 (Figure 3.5.2). Three 



data gaps are evident around CoRoT JD 3092, wiping out one transit completely. 

The cross-correlation function with a trapezoidal transit function is shown in Fig- 
ure 



3.5.2 The high stellar activity signal modulates the transit signals. 



The orbital period of the planet is 8.886593 ± 0.000004 days according to Fridlund 



et al. (2010). The rotational period was estimated with 6.4 ± 0.5 days differing slightly 



from the spectral analysis (6.9 ± 0.9 days). Looking closely on th autocorrelation func- 



tion (Figure 3.5.2) we clearly see that the transit-peak and the stellar rotation coincide 
after « 44.8 days. This is close to a 7 : 5-resonance between the stellar rotation and the 
orbital period. We might see here the evidence of an interaction between the planet 
and its host star. 
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Figure 3.14.: Cross-correlation function of CoRoT-6 with a trapezoidal function. 
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Figure 3.15.: Autocorrelation function of CoRoT-6. The orbital period of the planet 
is indicated by a dashed red line, the orbital period of the star by a 
dashed green line. The 7 : 5 resonance is marked with a blue dotted 
line. 
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Figure 3.16.: Lightcurve of CoRoT-2 with the transits of CoRoT-2b. The stellar 
activity signal is clearly visible and exceeds the depth of the transit at 
several times. 



3.5.3. Stellar Activity on CoRoT-2 

CoRoT-2b is a heavy and inflated planet orbiting an active G7V-star in a 1.7429-day 



orbit (Alonso et al. , 2008). Since its discovery, it has been an interesting object because 



of its unique features. The stellar activity signal is heavily modulated and changing 
over a period of ~ 30 days, which is a clear indication of spot evolution on this star 



(Figure 3.5.3). The magnetic activity of the star has been thoroughly investigated by 



Lanza et al. (2009) 



3.5.4. CoRoT-4 



We want to investigate CoRoT-4b where the planet's orbit of 9.20205 ± 0.00037 days 



and the stellar rotation 8.87 ± 1.12 days are closely related (Aigrain et al. , 2008). 



The spectral type of the host star is a G2V dwarf star which was determined from 
spectroscopic observations as well as the stellar radius i?* 



1.17±g;Bs (|Moutou et al. 



2008). The lightcurve can be seen in Figure 3.17 



The stellar rotation period was determined by calculating the autocorrelation-function 



see Aigrain et al. , 2008, for details). In order to analyze the stellar activity, the transits 
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Figure 3.17.: Lightcurve of CoRoT-4 with the transits of CoRoT-4b. If we compare 
to CoRoT-2, we see that the stellar noise is much higher. 



signatures have been taken out of the lightcure with the missing data being interpolated 



linearly. The autocorrelation function of the lightcurve is shown in Figure 3.19 



We can see the maxima of the stellar rotation signal at the determined period given 



by Aigrain et al. (2008), but we also see a signal, though less significant at half the 
period at ~ 4.5 days. This signal reaches its maximum at ~ 39 days, which can be 
explained with spot evolution on the rotating host star. 



By applying the transit detection algorithm described in Section |2.6| by using a 
trapezoid as a template for the cross-correlation function, we get a 6.51<r-detection 
for the transit signal. 

By applying a MCMC-method for the parameters duration and impact-factor of the 
trapezoid, the signal of the ccf was maximized resulting in an optimal duration of the 
template transit of 6.15 hours. 



From Equation |1 .13 and a total observation of 57.736296 days we can observe a 
maximum period P = 19.245432 days by assuming three fully observed transits. For 
a solar like star we would expect tx = 4.864 hours by applying Equation L9 An 
explanation for the difference between the maximum transit duration and the maximum 
of the ccf lies in the assumption for the solar like star, whereas CoRoT-4 was reported to 
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Figure 3.18.: Autocorrelation function (acf) of CoRoT-4b with transits (black) and 
without transits (red). The stellar rotation period is indicated by a 
vertical dashed line and its uncertainty with red dotted lines. The orbital 
period of the planet is marked with a black vertical dashed line. 
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Figure 3.19.: Crosscorrelation function (ccf) of CoRoT-4b. The ccf has been scaled 
with the standard-deviation of the ccf. 
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3.5. Exoplanets around Active Stars 
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Figure 3.20.: Dependency of the crosscorrelation function (ccf) in the impact factor 
(left) and the duration (right). The expected maximum duration of 
tx = 4.864 hours is marked by a vertical dashed line. 



have -R* ~ 1.17i? . We also have to keep in mind that the used template is a trapezoid 
and not a synthetic transit-shape, which would have introduced more free parameters 



(the limb darkening coefficients). As seen in Figure 3.20, the detection of the transit is 
insensitive to the impact factor, with a minor difference of 0.025a between the optimum 
b = 0.757 and the minimum. So the shape (triangular, trapezoidal or box-like) is less 
critical. 
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Stellar activity increases the difficulty to observe transiting planets and to measure their 
masses using radial velocity methods. Stellar spots change the size of the measured 



transit and increase the uncertainty by 3% (Czesla et al. , 2009). Besides that, the 
color temperature of the star is decreased by large spot groups which can lead to false 
characterization of the host stars properties and its mean density. 

The mean activity level of all spectral types is 377 ± 100 ppm, which is more than the 
average sun at 244.13 ppm. Even if the background noise of 144.3 ppm is of instrumental 
origin, the majority of stars is more active than the sun. Late type stars seem to be 
more active than earlier types. 

The radial velocity (RV) method is a proper tool to detect planets around main 
sequence stars ( Frink et al.| 2002). High stellar activity, however, limits the detection 
capabilities due to stellar granulation noise. This is a well known fact for radial velocity, 
but not for transiting systems. The problem of disentanglement has been discussed 



e.g. by Boisse et al. (2011). CoRoT-2b (Alonso et al. 2008) is one of the most famous 
examples for an extrasolar planet in orbit of an active star. CoRoT-6b is a good 
candidate for a planet-star interaction, showing a 7 : 5-resonance between the orbital 
cycle of the planet and the rotation period of the star. 



Planetary transits can be used to identify the stellar rotation period (e.g. Silva- 



Valio, 2008). Using the information from the transits, especially the inclination can 



be constrained, and the position of the spot on the star can be assumed. There are 



still drawbacks, as described by Silva-Valio (2008), since a spot can not be clearly 
identified and there is an ambiguity in the observation of HD209458, where the period 
could either be 9.9 or 11.4 days. Nevertheless this method is a valuable alternative to 
spectroscopically measured line broadening. 



As seen in Figure |3.20[ the detection of the transit is independent of the impact 
factor, regardless of the stellar activity. So the shape (triangular, trapezoidal or box- 
like) does not favour the detection. This fact is used the the transit detection algorithm 
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by Kovacs et al. (2002) (see Section 2.2). 

A robust filtering algorithm is used to clean the lightcurve from instrumental effects 
and to suppress the stellar activity signal. Wiener filtering is not suitable for filtering 
stellar activity since it fails to estimate the noise background. Usually the subtraction 
of a median average or application of a Butterworth filter is sufficient. In rare occasions 
where pulsations or oscillations pollute the lightcurve, pre-whitening is most effective 
to enhance the transit detection. 



The application of the filter by Savitzky and Golay ( 1964 ) is essential to the quality of 
the data, but has a severe impact on the transit-shape transforming it into a Gaussian. 
Alternatives have to be considered to ensure a robust outlier-detection and removal. 
Generally simple algorithms like e.g. a three-point median filter is not sufficient. 

Nevertheless, a well designed filter system can enhance the transit detection and 



suppress the stellar activity signal (Weingrill et al. 2010). 

Stellar variability can be easily identified by investigating the lightcurve using auto- 
correlation. Using radial velocity measurements additionally helps to constrain the 
rotational period of the star. 

Transit signals have to be expected in the domain of stellar rotation periods. CoRoT- 
2b and CoRoT-6b are strikingly examples for these problems. Nevertheless, the transit 
signal usually has sufficient power at higher harmonics to be identified. 

Except for multi-band observations like the three color channels of CoRoT stellar 
flares are almost impossible to be distinguished from instrumental effects like cosmic 
ray hits. Amplitudes and decay-rates and post-event-levels might give indications to 
detect a flare. Nevertheless such rare events do not harm the transit detection. 

A magnetic or gravitational interaction between a planet and its host star is most 
likely as demonstrated on CoRoT-6b. Tidally locked planets or planets near the Roche- 
Limit are candidates for direct or resonant interaction. 

Stellar pulsations like p- modes do not impair the detection of a transiting planet. 
This might be the case for radial pulsators like e. g. Cepheids. We can overcome this 
problem either by modelling pulsations or by using additional constraints like the color 
information. The same holds for binary stars, which were not discussed within the 
frame of this work, where planetary transit signals can be detected after removing the 
binary transit signals. 



We are now able to extend the values given in Table 1.2 with semi- major axis and 



their respective periods for the habitable zone taken from Kasting et al. (1993) with 
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the transit depth estimated for a Super-Earth (Table 4.1). 



Table 4.1.: Spectral types of dwarf stars taken from Cox (2000) with their respective 
masses M* in solar masses and solar radii . The table has been extended 
with the expected semi-major axis a and the Keplerian period P as well as 
the expected transit depth AF in ppm for a transiting Super-Earth with 
two Earth-radii and the approximate transit duration tx in hours. 



Spectral Type 


M* 


-R* 




a 




P 


AF 


t T 


A5V 


2.0 


1.7 


2 


81 


3 


34 


116 


52.34 


FOV 


1.6 


1.5 


2 


20 


2 


58 


149 


45.68 


F5V 


1.4 


1.3 


1 


82 


2 


08 


199 


38.48 


GOV 


1.05 


1.1 


1 


13 


1 


17 


278 


29.60 


G2V 


1.0 


1.0 


1 


03 


1 


05 


336 


26.39 


G5V 


0.92 


0.92 
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87 


397 


23.43 


KOV 


0.79 


0.85 





66 





61 


465 


20.20 


K5V 


0.67 


0.72 





48 





41 


648 


15.80 


MOV 


0.51 


0.60 





28 





21 


933 


11.51 


M2V 


0.40 


0.50 





18 





12 


1344 


8.62 


M5V 


0.21 


0.27 





08 





05 


4608 


4.37 



Given the level of stellar activity determined from CoRoT the stellar activity of early 
types might prevent the successful detection of a Super-Earth in the habitable zone. 
The sun as a G2V star within an average activity level is the limiting case. 



In rare occasions flares are observed during transits (Bentley et al. , 2009). But 
photometric methods that measure the integral flux of the star can not discern between 
a temporal phenomenon like a flare on the stellar disk or a spatial phenomenon like a 
stellar spot. 

Future missions will utilise new methods like transit spectroscopy or observation in 
different wavelength regimes to identify stellar activity and its interaction with the 
hosted planet. 
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2.11. Power spectrum of CoRoT-6 with signal and noise levels estimated with 
a polynomial. The signal is indicated with a dashed line, the noise is 
indicated by a dotted line. Both signal plus noise are represented by the 
dashed dotted line. We see that the signal is underestimated, while the 
level of noise is somewhat accurate EE] 



3.1. 



3.2. 



3.4. 



3.5. 



3.6. 



3.7. 



3.8. 



High frequency part of the power spectrum of a F9V star CoRoT LRaOl E2 0205. 
The orbital frequency of 0.323 mHz and its harmonics are indicated by 
dashed vertical lines 

A simulated transit (black solid) and the results of the filtering process. 
The red crosses indicate the 512 second oversampled data, green after a 
median filter was applied and blue after applying |Savitzky and Golay 



(1964) 



3.3. First white-light flare observed with a TSI (red diamonds) instrument. 



The NOAA/ GOES X-ray flux is shown for comparison (blue squares) . ( Kopp 



et al, 2004 Woods and Kopp, 2005) 



Total Solar Irradiance measured with the TIM/SORCE instrument with 
a 6 hour cadence. The data covers a little less than a solar cycle. High 
activity signals correspond to a higher output of the sun. The largest 
event in 2003 was caused by a solar spot 

Correlation between the Total Solar Irradiance measured with the TIM/SORCE 
instrument and its variance over a timescale of 150 days (blue) and over 
30 days (red). The lengths have been chosen in accordance with the 
CoRoT observational run lengths 



E3 



Power spectrum of a F9V star CoRoT LRaOl E2 0205. The spectral 
power decreases from low frequencies to high frequencies. A small peak 
at the one day period is visible, as well as stronger peaks for the orbital 
period at 0.323 mHz and its harmonics (see also Ballot et al. , 2011). 



The granulation noise is visible between 0.1 and lmHz. The p-mode 
oscillations are centred around 4 mHz 

Stellar noise levels within 6 hours on the CoRoT run LRaOl. The solar 
noise level is shown in red indicating the minimum and maximum noise 
signal during the solar cycle 

Periods in cycles-per-day (cpd) versus amplitude of all stars in the CoRoT 
run LRaOl. The red line corresponds to a transit depth of Jupiter, the 
green line corresponds to a Neptune and the blue line to an Earth-sized 
transit depth 



EH 



ED 
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3.9. Partial power spectrum of a F9V star CoRoT LRaOl E2 0205. The 
spectral power decreases from low frequencies to high frequencies. A 
small peak at the one day period is visible, as well as stronger peaks 



for the orbital period at 0.323 mHz and its harmonics (see also Ballot 



et al. 2011). The granulation noise is visible between 0.1 and lmHz. 



The p-mode oscillations are centered around 4 mHz [72] 

3.10. Solar image form 6 June 2011 taken with the SDO/AIA instrument at a 
wavelength of 450 nm. Four different quadrants are visible in the image, 
which might be caused by instrumental effects [73] 

3.11. Simulation of a transiting planet with two Earth radii. The black sold 
line is the simulated transit measurement. The red line is the fit of a 
model by Mandel and Agol (2002). The solid green line are the residuals 

of the fit [7J 

3.12. Enlarged residuals from Figure 3.11 The diameter of the planet is rep- 
resented of the vertical dashed lines. The center of the star is indicated 

by the vertical dashed dotted line [75] 

3.13. Lightcurve of CoRoT-6 with the transits of CoRoT-6b. In the last 
quarter the stellar activity signal is in the order of the transit depth. . [75] 

3.14. Cross-correlation function of CoRoT-6 with a trapezoidal function. . . . [77J 

3.15. Autocorrelation function of CoRoT-6. The orbital period of the planet is 
indicated by a dashed red line, the orbital period of the star by a dashed 
green line. The 7 : 5 resonance is marked with a blue dotted line. . . . [77] 

3.16. Lightcurve of CoRoT-2 with the transits of CoRoT-2b. The stellar activ- 
ity signal is clearly visible and exceeds the depth of the transit at several 
times [78] 

3.17. Lightcurve of CoRoT-4 with the transits of CoRoT-4b. If we compare to 
CoRoT-2, we see that the stellar noise is much higher [79] 

3.18. Autocorrelation function (acf) of CoRoT-4b with transits (black) and 
without transits (red). The stellar rotation period is indicated by a 
vertical dashed line and its uncertainty with red dotted lines. The orbital 
period of the planet is marked with a black vertical dashed line [8U] 

3.19. Crosscorrelation function (ccf) of CoRoT-4b. The ccf has been scaled 
with the standard-deviation of the ccf [8U] 

3.20. Dependency of the crosscorrelation function (ccf) in the impact factor 
(left) and the duration (right). The expected maximum duration of tx = 
4.864 hours is marked by a vertical dashed line EH 
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A.l. The drizzling algorithm 

i function drizzle, t, y, err=err, LOCATIONS^! o c a t i o n s , B I N S I Z E=B I N S I Z E , $ 

_EXTRA = extra 

; I a st bin is empty ! 

if N_ELEMENTS(t) LE 1 then $ 

message , ' insufficient u times 
e if N_ELEMENTS(t) NE N_ELEMENTS( y ) then $ 

message , ' t i mes u and u d a ta u not u of u same u I e n gt h 
IF BINSIZE EQ 0.0 then begin 

message , ' binsize u zero ! ' , / continue 

err = REPLICATE(STDDEV( y ) , N_ELEMENTS( y ) ) 
ii locations — t 

retu rn , y 
e n d i f 

if not Keyword_set(NBINS) then $ 

NBINS = CEIL((MAX(t ,/NAN)-MIN(t , /NAN))/ BINSIZE) 
i6 h = histogram (t , REVERSEJNDICES = r , / NAN , LOCATIONS = locations , $ 
BINSIZE = BINSIZE, _EXTRA = extra) 
result = DBLARR( NBINS, /NOZERO) 
err = DBLARR( NBINS , /NOZERO) 
result [*] = I VALUES . F_NAN 
21 for i = 0L, nbins— 1L do begin 

IF R[ i ] NE R[ i+lL] THEN begin 
ybin = y [R[R[ I ] : R[ r+lL] 1L ] ] 
result [ i ] = MEAN(ybin ,/NAN) 
; WARNING: stddev ignores the nan keyword I 
26 err[i] = STDDEV( y bi n , /NAN) 

; alternatively: 

; err [ i ] = RMS( ybin ) 

e n d i f 
e n d f o r 

3i if N_ELEMENTS( result ) NE N_ELEMENTS( locations ) then $ 

message , ' Warn ing : u u n eq u a I u a r r a y s u i n u r e s u 1 1 ' , /CONTINUE 

retu rn , result 
end 

A. 2. Cosine Fit 
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PRO COSFUNC, X, A, F, pder 

A[0] ... xl Phase in days 

1 ... x2 duration in days 

2 ... x3 amplitude 

3 . . . x4 offset 
F = cos((X+ A[0])*2.0*! DPI/A [ 1 ] ) * A[2]+A [3] 
IF N_P ARAMS ( ) GE 4 THEN $ 

pder = [ [-2.0*!DPI*A[2]* sin (2 . 0* ! DPI *(X+A [0] ) / A[ 1] ) / A [ 1] ] , $ 

[2*! DPI*A[2]*(X + A[0])* sin (2 . 0* ! DPI *(X+A[0] ) /A [1] ) / (A [1] ~2 )] , $ 
[cos (2.0*! DPI*(X+A[0])/A[1])] , $ 
[ replicate (1.0 , N_ELEM ENTS ( X ) ) ] ] 
RETURN, RMS(I1 - (cos((tl+ P0 [ 0] ) * 2 . * ! DPI/P0 [ 1 ] ) * P0 [2] + P0 [ 3 ] ) ) 
;D[cos((t + xl)*2.0*Pi/x2)*x3+x4,xl] 
d/dxl = -2.0*! DPI*x3*sin (2.0*! DPI *( t+xl )/x2 )/ x2 
d/dx2 = 2*!DPI*x3*(t + xl )* s i n ( 2 . * ! DPI *( t+xl )/ x2 )/ ( x2^2) 
d/dx3 = cos(2.0*! DPI *( t+xl )/x2 ) 
d/dx4 = 1.0 
END 

20 FUNCTION COSFIT, X, Y, Weights, A, Sigma 

RETURN, CURVEFIT(X, Y, Weights, A, Sigma, FUNCTION_NAME= ' COSFUNC ' ) 
END 

A. 3. Fourier Autocorrelation 



FUNCTION ACORR, x, NORMALIZE=NORMALIZE 

f 

NAME: 

ACORR 

PURPOSE: 

Compute the fourier autocorrelation function 
CALLING SEQUENCE: 
ACORR, x 

INPUTS: 

x data to a u t oco r re I a t e 

OPTIONAL INPUTS: 
none 

OUTPUTS : 



result 
KEYWORD PARAMETERS 
NORMALIZE 



a u t oco r re I a te d data 



performs normalization in respect of the finite length 
of the dataset. 
MODIFICATION HISTORY: 

Version 1.0, 04.03.2011, Joerg Weingrill SRI/AAS. 

(joerg . weingrillOoeaw . ac . at) 

n = N_ELEMENTS( x ) 
23 k = DINDGEN ( n )/ n 

xl = [x , REPLICATE(MEAN(x) , N_ELEM ENTS ( x ) ) ] 
fftx = fft (xl,-l, /DOUBLE) 
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ac = REAL_PART( fft ( fftx*conj ( fftx ) ,1)) 

if KEYWORD_SET( NORMALIZE) then RETURN, ( ac/ac [0] ) / ( 1 . - k) else $ 
RETURN, (ac/ac[0]) 

END 

A. 4. Fourier Cross-correlation 



FUNCTION CCORR, x,y 

; calculate the c ross — c o r r e I a t i o n between x and y 
RETURN, REAL_PART( fft (conj ( fft (x, -1))* fft (y , -1) ,1)) 
4 END 

A. 5. Fourier Interpolation 



; compute the fourier interpolation of a series 
FUNCTION FFTFUNCT, P, X=X, Y=Y 
j = WHERE(~ FIN ITE ( Y) ) 
m = N_ELEMENTS(P) 
5 Yl = Y 
Yl[j] = P 

RETURN, REPLICATE(TOTAL(ABS( FFT( Yl ) ) ) ,m)/m 
END 



function FINTERPOL, V, X, XOUT 

if N_ELEMENTS(XOUT) EQ then begin 

; if we do not have XOUT, assume non finite values to be interpolated . 
i = WHERE(~FIN ITE (V) ) 
is XOUT = X[ i ] 
e n d i f 

; use linear interpolation as a first guess 
; usually zeros would do the same 
P0 = INTERPOL(V.X.XOUT) 
20 FA = {X:x, Y:v} 

; MPFIT is far more robust then AMOEBA algorithm 
R = MPFIT( 'FFTFUNCT' , P0 , f u n ct a rgs=FA , FTOL=le -5,MAXITER=10) 
RETURN, R 
end 

A. 6. Trapeozid Function 

i function trapezoid , x, xO , height , duration , impactfactor 

; + 



/ \ 

/ \ 

12 3 4 

xO ... time of mid trapez (t3— 12)/2 
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height ... hight of trapez 
duration ... t4-tl 
impactfactor ... ( t3-t2 )/( t4-tl ) 
flatpart = (t3— 12) = d u ra t i o n * i m p a c tf a ct or 

n = N_ELEMENTS( x ) 
; constraints 

if impactfactor GE 1.0 then impactfactor = 0.999 
i6 if impactfactor LT 0.0 then impactfactor = 0.0 
duration = ABS( d u ration ) 

y = DBLARR(n) 

if (height eq 0.0) then return ,y 
2i dt = (x[n-l]-x[0])/(n-l) 

duration2 = d u ra t i o n / ( 2 . d ) 
xl = DBLARR(6) 
yl = DBLARR(6) 

26 xOdt = xO ; /dt 

tl = (xO - duration2 ) 

t2 = (xO — ( d u ra t i o n 2 * i m p a ct f a cto r ) ) 

t3 = (xO + ( d u ra t i o n 2 * i m p a ct f a cto r ) ) 

t4 = (xO + duration2 ) 
31 xl = [tl -1.0 , tl , t2 , t3 , t4 , t4 + 1.0] 

yl = [0.0 , 0.0, height, height, 0.0, 0.0] 

y = INTERPOL( yl , xl , x) 

retu rn , y 
end 



A. 7. Trapeozid Fit 



FUNCTION TRAPEZFUNC.X , P 

return, trapezoid(x, P[0], P[l], P[2], P[3]) 
END 

pro mpfindtrapez , x, y, err , x0=x0 , height=height , d u ration=d u ration , $ 

impactfactor=impactfactor , BESTNORM=BESTNORM , PERROR=PERROR 
P0 = [xO , height , duration , impactfactor] 

R= mpfitfun ( 'TRAPEZFUNC , x, y, err, pO , $ 

PERROR=PERROR , BESTNORM=BESTNORM , /NaN, /QUIET) ; 
if n_elements ( r ) gt 1 then begin 
xO = R[0] 
height=R[l] 
duration =R [ 2 ] 
i m p a ctf a c to r=R [ 3 ] 
e n d i f 
end 
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A. 8. Root Mean Square 

FUNCTION RMS, X 

; compute the root mean square of the vector X 
3 RETURN, SQRT(MEAN(X^2 ,/NAN) ) 
END 

A. 9. Folding of Data 

ipro folded.t, I, period, t f=tf , I f =1 f , ferr = ferr 
tO = t [0] 

tf = (t-tO) mod Period 
i = SORT( tf ) 
tp = tf [i] 
6 lp = I [i] 

k = max(t-t0 ) / Period 

If = drizzle(tp, lp , err=ferr , loc=tf , Bl NSIZE=med ia n ( d e r i v ( t ) ) ) 
end 

A. 10. Period Testing 

i FUNCTION PERIODFIND.P0 

COMMON PERIODFUNC, ct , cl 

; fold mod P0 
tf = ct - ct [0] 
P0[0] = P0[0] < max(tf)/3.0 
e P0[0] = P0[0] > 1.0/max(tf) 
tf = tf mod P0[0] 
i = SORT( tf ) 
tp = tf [i] 
lp = cl [i] 

ii n = CEIL(0.25*P0[0]*86400.0d/512.0d) 
if n GT 1 then begin 

;tpl = P0[0]*DINDGEN(n)/(n-l) 
; drizzle with new period 
dl = drizzle(tp, lp, e r r=s , I o c=tp 1 ) 
i6 ; return rms of stddev in drizzle bins 

r = max(s) 
end if else r = RMS( cl ) 
RETURN , r 
END 

21 

pro periodtest 

outpath = 'c:\Users\jwe. SATGEO\ Documents \Dissertation\thesis\figures\' 
loadlightcurve , 'Q:\SRa02\data\SRa02_El_0359. fits ' ,t , I ,$ 
fluxdev=ld , name=name ; Cepheid 
26 n = N_ELEMENTS( I ) 
t = REFORM ( t [ 1 ,*]) 
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lr = REFORM ( I [0 ,*]) 

[g = REFORM ( I [1 ,*]) 

lb = REFORM ( I [2 ,*]) 
si Ird = REFORM ( Id [0 ,*]) 

n = (max(t)-min(t))*86400./512. 

tl = FINDGEN ( n ) * ( max( t) — min ( t ) ) / n + min(t) 

Irl = drizzle(t,lr ,nbins = n,lo c=t 1 ) 

Igl = d r i z z I e ( t , Ig , n b i n s = n,loc=tl) 
36 Ibl = drizzle(t,lb,nbins = n,lo c=t 1 ) 

Irdl = d r i z z I e ( t , I rd ,nbins = n,loc=tl) 

i = WHERE( I r 1 GE . 1 ) 

Irl = Irl [i ] 

Irdl = Irdl [ i ] 
4i i = WHERE( Igl GE 0.01) 

Igl = lgl[i] 

i = WHERE( Ibl GE . 1 , n ) 
Ibl = Ibl [ i ] 
tl = tl[i] 

46 mr = -2. 5*ALOG10( Irl )*0. 937906 + 25.5486 
mg = -2. 5*ALOG10( Igl )*0. 703167 + 21.9155 
mb = -2. 5*ALOG10( Ibl )*0. 963139 + 25.6203 



P = 8 ;t2[i [j]] 
si II = Irl + Igl + Ibl 

II /= MEAN( II ) 

COMMON PERIODFUNC, ct , cl 

ct = tl 

cl = II 
56 n = 500 

v = dblarr(n) 

p = dblarr(n) 

for i=0,n— 1 do begin 

;p[ i ] = 6.5 + i *(1.0/n) 
ei p[ i ] = 0.5 + 10.* i *(1.0/n) 

v[i] = pe r i o d f i n d ( p [ i ] ) 

e n d f o r 



mv = min ( v , i ) 
ee p r i n t , p [ i ] , mv 

PS_Start , FILENAME=outpath+' periodtest.eps', /European, FONT= 1 , TT_FONT= ' A r i a I ' , 
LANDSCAPES, INCHES=0, XSIZE = 15, YSIZE = 10, XOffset=0, YOffset=0, /ENCAPSULATED 
plot ,p,v,/ xs ,/ ys , title =name, xt i 1 1 e = ' period u [days] ' , $ 
ytitle='max(sigma) ' ,/nodata , yrange = [0. 0,0.1] 
7i oplot , p , v , co I o r=FSC_COLOR( ' red ' ) 



PS_END 

P = p[ i ] 

1 1 — = mean (II) 

76 ;findperiod,tl,ll, p=p0 , c h i sq=c h i sq 
n = N_ELEMENTS( II) 
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fdr = FFT([REPLICATE(0.0d,n/2),[l , REPLICATE ( . d , n/2)] ) 
n2 = N_ELEMENTS( 1 1 ) 
ac = acorr(ll) 
si ac2 = A_CORRELATE( 1 1 , INDGEN(n)) 

freqs = 86400 . Od *FINDGEN ( n2/2 + l)/( n2 * 5 1 2 . Od ) 
periods = 1.0/freqs 

plot .periods ,ABS(fdr [ : n - 1] )~2 ,/xs , / ys , / y log , $ 
xtitle = 'period u [days] ' ,xrange=[0.5,11.5] 
86 plot , tl [0 : n-l]-tl [0] , ac 

; print ,p0,chisq 
end 

A. 11. Prewhitening 

PRO gfunct , X, A, F, pder 
2 n = n_elements (x) 

F = A[0] * SIN(A[1]*2.*I DPI*X/n + A[2]) 
daO = SIN(A[1]*2.*! DPI*X/n + A[2]) 

dal = A[0] * COS(A[l]*2.*! DPI*X/n+ A[2])*2*l DPI*X/n 
da2 = A[0] * COS(A[l]*2.*l DPI*X/n + A[2]) 
7 IF N_PARAMS() GE 4 THEN $ 

pder = [[daO] , [dal] , [da2]] 

END 



12 FUNCTION prewit,x,y,maxiter, freqs = freqs, amps = amps, phas = phas 

; + 

;prewitning of function y(x) with maxiter iterations 

freqs = f 1 1 a r r ( m axi te r ) 
17 amps = f 1 1 a r r ( m axi te r ) 
phas = f 1 1 a r r ( m axi te r ) 

;freqs_sig = f 1 1 a r r ( m axi te r ) 
;amps_sig = f 1 1 a r r ( m axi te r ) 
;phas_sig = f 1 1 a r r ( m axi ter ) 
22 n = N_ELEMENTS(y) 
dt = x[l]-x[0] 



for iter = 0, maxiter— 1 do begin 
27 f = (abs(fft(y)))[0:n/2-l] 

pha = (ATAN( fft (y) , /PHASE) ) [0: n/2-1] 
freq = f i n d ge n ( n / 2) / ( n * dt ) 

weights = re p I i c a te ( 1 . , n ) 

32 

mf = max( f , i ) 

mfi = n /(l./ freq [ i ]) 

A= [mf*2. , mfi , pha [ i ] + ((! DPU0.5) MOD I DPI)] 
status = 
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37 sigma = 

; plot ,x,y, psym=3 

yfit = CURVEFIT(X, Y, weights, A, SIGMA, FUNCTION_NAME= ' gfunct ' , status=status ) 

; gf u net , x , a , yfit 
;oplot,x,yfit , color =1 
42 amps [iter] = a [0] 

freqs [iter] = a [1] 
phas[iter] = a[2] 

if status eq then y — = yfit else message , 'no u fit I 
; print , a, sigma 
47 endfor 

retu rn , y 
END 
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